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SUMMARY 

The various steps that could control the vaporization rate of a material 
are discussed. These steps Include the actual vaporization, flow rate of 
matrix gas, chemical reaction, gas diffusion, and solid state diffusion. The 
applicable equations have been collected from diverse appropriate sources, and 
their use Is explained. Rate equations are derived for conditions where more 
than one step Is rate controlling. Calculations are made for two model 
materials: rhenium which vaporizes congruently, and tantalum carbide which 

vaporizes Incongruently. The case of vaporization under thermal gradient 
conditions Is also treated. The existence of a thermal gradient In the 
reslstojet means that the vaporization rate of a material may be only one 
thousandth of that predicted under Isothermal conditions. Calculations show 
that rhenium might have a 100 000 hr lifetime at temperature In a 2500° C 
reslstojet. Tantalum carbide would have a life of only 660 sec under similar 
conditions . 


INTRODUCTION 

Over the past years, there has been frequent need to calculate the 
expected vaporization (sublimation) rates for materials to be used In various 
design concepts. Recently, we have made such calculations for a 2500° C ver- 
sion of the reslstojet. The theoretical equations are well established, but 
their application to specific conditions Is sometimes a problem. It occurred 
to us that It would be helpful for future use to assemble Into a single docu- 
ment (this report) these equations as well as effective ways to use them In a 
practical vaporization problem. Because of the large amount of Information to 
be covered In this report, we have divided It Into four parts; a preview of 
each part follows: Parts I and II are concerned with the understanding of the 

principles and derivation of rate equations. Parts III and IV are oriented 
toward practical problems. 

In Part I, we first discuss the general vaporization concepts and define 
several ways of expressing vaporization rates. In Part I, we also present 
equations for calculating vaporization rates If a single kinetic step controls 
the process. The simplest calculation of vaporization rates, of course. 
Involves the (Langmuir) equation which gives the rate for the vaporization 
step. Itself. This equation Is applicable when a material Is exposed directly 
to vacuum. In actual practice, this normal vaporization (nv) step may not be 
rate controlling. When vaporization occurs within a partial enclosure, the 
effective rate Is slowed down, and the vapor species builds up In the con- 
tainer; the rate becomes controlled by the rate of escape from the container 
which Is determined by the area of the escape hole. When a matrix gas Is 
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present In the system, there are several steps that might be rate control- 
ling: the flow of the matrix gas, a chemical reaction, or the diffusion of 

the vaporizing species through the matrix gas. Moreover, If the material 
being evaluated vaporizes Incongruently, solid state diffusion may become rate 
controlling. Incongruent vaporization refers to materials for which the ratio 
of the rates of vaporization of the various components of the material Is not 
equal to their ratio In the solid material. For such materials the composi- 
tion of the solid material changes as vaporization progresses. Note that this 
report Is not concerned with vaporization Induced by laser, electron, or Ion 
bombardment. 

The experimental conditions determine which step Is rate controlling. 
Aprlorl, one cannot know which steps In a particular overall process are 
Important. Therefore, In order to discuss the vaporization process, we had to 
direct the discussion and calculations toward a specific system. We chose the 
reslstojet; It Is a good model system because most of the aforementioned 
potentially rate-controlling steps might apply. As model materials we chose 
rhenium (Re) and tantalum carbide (TaC). These are good choices because the 
vaporization behavior of the one contrasts the other: the Re vaporizes con- 

gruently, and TaC vaporizes Incongruently. Moreover, chemically TaC Is re- 
active while Re Is relatively unreactive under the conditions of Interest. 

The principles applied to Re and TaC In this report apply equally well to 
ceramics, polymers, alloys, composites, and all other materials. 

We begin Part II by giving the equations that show how kinetic steps can 
be combined to give kinetic rates that are controlled by two or more Individ- 
ual steps. Via these equations, we have derived the equations which give the 
vaporization rate when the steps discussed In Part I are coupled. Some of the 
equations are relatively complex, being dependent on both time and location In 
a system. Calculation of vaporization of both Re and TaC were again made; 
dependence of rates on time and location are Illustrated graphically. 

It should be noted that the actual rate of vaporization In a complex 
physical system can only be obtained experimentally. A calculated value Is 
generally an equilibrium or steady state value which may never actually be 
reached; the calculated value should be considered to be merely a maximum. 

Also In this report, several approximations are used, both In parameters (such 
as vapor pressure) and In the derivations of the equations. Moreover, Ideal 
behavior Is assumed In most cases so that the equations will be more widely 
applicable. Some of the calculations are based on the parameters specific to 
our concept of the reslstojet, so that the final conclusions as to the nature 
of the rate-controlling step may not be applicable to reslstojet concepts with 
quite different parameters. 

In Part III, we analyze the kinetic steps as to which are most Important 
in an Isothermal system. Our reasoning Is then applied to the reslstojet con- 
cept, assuming Isothermal conditions. Calculated vaporization rates for Re 
and TaC are again given. Other topics discussed In Part II are the effect of 
heating and cooling periods and the use of a sacrificial vaporizing material. 

In the beginning of Part IV, we discuss the vaporization rate as a func- 
tion of temperature. Using two different approaches, we derive rate equations 


2 



for the condition of a temperature gradient. When appropriate, the vaporiza- 
tion rates for Re and TaC were calculated as a function of gradient, location, 
and time. The significance of a temperature gradient In the reslstojet Is dis- 
cussed. Finally, the expected lifetimes of Re and TaC are calculated using 
all the pertinent kinetic steps and assuming the existence of a 10° C/cm gra- 
dient. From a vaporization viewpoint. Re appears to be a satisfactory material 
with a potential lifetime of 100 000 hr while TaC does not. 

Our findings are summarized at the end of this report In the Synopsis. 

As an aid to the reader In the use of the equations presented, we have Included 
a Glossary of symbols In the Appendix. In order for this report to be useful 
to people of all fields for designs In which vaporization may be a problem, 
simple language Is used and concepts are repeated where necessary. Of course, 
all sections cannot be of uniform simplicity. The source and use of equations 
are explained, but detailed derivations of equations are minimized. When der- 
ivations are given, they are Intended as a method of explaining the use and 
meaning of the desired equation. 


THE RESIST0JET 

At this point It Is worthwhile to give a little discussion about the re- 
slstojet (ref. 1). The parameters of one design are given because some of the 
vaporization equations In this report require the use of such values. 

The original concept of the reslstojet Involved merely the resistive heat- 
ing of a gas to be expelled from a spacecraft to get added thrust. A more 
recent concept Includes an additional Idea, that of converting waste material 
Into a gas and then expelling It from the reslstojet to get extra thrust. 

Thus, the newer concept of the reslstojet would solve two problems. The design 
for the reslstojet Is continually changing. The schematic given In figure 1 Is 
sufficient for our discussion of the vaporization problem. 

The best gas for the operation of the reslstojet In terms of thrust Is 
hydrogen ( H 2 ) - However, other gases that may be expelled are: water vapor, 

carbon dioxide, carbon monoxide, nitrogen, ammonia, and hydrazine; these gases 
give less thrust than hydrogen and. In general, are more corrosive. The tem- 
perature proposed for the reslstojet Is 2500° C (2773 K). At present, the 
desired time for operation of the reslstojet at this temperature Is 1000 hr; 
but the desired time will probably Increase. In our discussion of vaporiza- 
tion, we are Interested In the hottest portion of the system, the place where 
the vaporization rate Is expected to be the greatest. This Is the final tube 
before the nozzle (fig. 1). This hottest portion of the system Is also the 
place that Is most susceptible to chemical reaction with the flowing gases. 

This tube might possibly be made from either rhenium or tantalum carbide be- 
cause both of these materials have high melting points and are electrical con- 
ductors which enables direct electrical heating. 

The calculations In this report are based on the following parameters for 
the reslstojet. The tube's length (L) Is assumed to be 10 cm; Its Inside diam- 
eter (2R‘), 1 mm; and Its wall thickness (S), 0.1 mm. A value of 3.16 atm Is 
taken as the pressure ( P(m) ) of the flowing matrix gas. The pressure drop 
along the tube Is presumed to be negligible. The flow rate ( F ' ) Is 5000 stand- 
ard cmVmln (268 cmfysec at 2773 K and 3.16 atm pressure). 
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PART I: SINGLE STEP RATE-CONTROLLING PROGRESS 

General Concepts 

In the actual vaporization process, the molecules (or atoms) do not just 
fly off the crystal lattice of the solid. Some of the molecules on the surface 
are undoubtedly mobile, not bound to a particular position of the crystal lat- 
tice but free to move over the entire surface of the material. If molecules 
on the surface were not mobile, recrystalllzatlon would be Impossible. The 
mobile molecules that have sufficient energy can leave the surface through 
vaporization. The number which vaporize per unit should be proportional to the 
activity (number) of the mobile molecules on the surface for any given temper- 
ature. This Is a physical picture of vaporization but does not lead to an 
equation for the calculation of the rate of vaporization. 

To obtain the desired equation, one must start with a picture of the mol- 
ecules In the vapor above the solid and assume they are In equilibrium with 
the molecules In the solid. Equilibrium means that the number of molecules 
vaporizing In a given time Is equal to the number condensing (returning to the 
solid from the gaseous state) In the same period of time. Thus, the number of 
molecules In the solid does not essentially change at equilibrium. The rate 
of condensation, N', Is equal to the rate at which molecules strike a unit 
surface area In a unit time multiplied by the fraction, Q ' , that adheres. The 
resultant equation Is known as the Langmuir equation (ref. 2, p. 278). 

N' = 1.013X10 6 P Q'/sqr(2irRTM) = a P/sqr(M) (1) 

where N 1 Is the rate In terms of moles/cm 2 sec, a mole being the number of 
molecules that has a weight In grams equal to the molecular weight. The P 
Is the pressure In atm; the factor, 1.013x10^, converts the pressure to 
dynes/cm 2 ; H Is the molecular weight of the species In grams; R Is the gas 
constant (8.31 xlO 7 ergs/K); T Is the temperature In K; and the "a" Is equal 
to 1 .01 3x1 0^/sqr( 2irRT) which Is 0.842 g ^ 2 moles^ 2 cm 2 sec - ^ atm - ^ 
at T = 2773 K. 

The Q' In equation (1) Is called the sticking factor or accommodation 
coefficient and can have any value from zero to one. When Q' Is set equal 
to one, we again have the equation that gives the number of moles of a material 
striking a unit area In a unit time. It Is known as the Knudsen equation. The 
Knudsen equation Is derivable from the kinetic theory of gases and has been 
shown to be correct experimentally under near vacuum conditions. In this re- 
port, Q' Is always assigned the value of one. As alluded to before, under 
equilibrium conditions, the N' In equation (1) becomes the number of moles 
vaporizing from a square cm In a sec, P becomes the vapor pressure (P(v)), and 
M Is the molecular weight of the vaporized species. Thus, we now have an 
equation that Is applicable to the vaporization process under vacuum condi- 
tions. Under nonvacuum conditions, the net rate of vaporization Is determined 
by the difference between the vapor pressure and the pressure of the vaporized 
species In the surrounding space. Because this pressure can be a function of 
the distance (Y) along the enclosure (tube), we refer to It as P(Y). It fol- 
lows that 


N' = a (P(v) - P(Y))/sqr(M) 


( 2 ) 
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One can more easily comprehend the magnitude of the vaporization rate If 
It Is expressed In units of g/cm 2 sec. We will call this quantity, W, and 
the equation by which It Is calculatable Is obtained by multiplying equation 
(2) by the molecular weight of the vaporizing species. 

W = a (P(v) - P( Y) ) sqr(H) (3) 

When vaporization occurs Into vacuum, we will term the process normal vaporiza- 

tion and equation (3) becomes 

W(nv) = a P(v) sqr(M) (4) 

The total rate of vaporization can always be obtained by multiplying the aver- 
age rate per unit area by the Inside area of the tube (A(s) = 2irR'L with R' 

being the radius and L, the length). 

The recession rate, G, Is another way of expressing the rate of vaporiza- 
tion. It Is the rate at which the material decreases In thickness (cm/sec). 

Its value can generally be obtained by dividing W by the density of the mate- 
rial, D. Thus, 


G = a (P(v) - P( Y) ) sqr(M)/D (5) 

By using the room temperature density, one can compare G with the room tem- 
perature thickness. 

Throughout this report, we have endeavored to use a uniform set of sym- 
bols. However, because W always refers to the rate of vaporization, this Is 
equivalent to continually redefining It In terms of the rate-controlling step 
presently being considered. In a similar manner, N‘, and G are likewise 
redefined. Doing this emphasizes the concept of loss by vaporization as well 
as minimizes the use of subscripts and foreign alphabets. Another symbol that 
we use In this way Is Z, the ratio of the rate of vaporization under a partic- 
ular set of conditions (pc) to the rate under normal vaporization condition 
(nv). 


Z = W(pc)/W(nv) = G(pc)/G(nv) 
= N'(pc)/N'(nv) 


( 6 ) 


The final form of the equation Is only true If no chemical reaction occurs. 

The reason for this will become evident In a subsequent section; It Involves 
the difference between the molecular weight of the normally vaporizing species 
and that of the chemically formed species. The calculation of Z for each 
rate-determining step eliminates the necessity of always calculating three 
different values: W, G, and N'. Rates that are time or space dependent can- 

not always be expressed In terms of Z. Therefore, at appropriate places In 
this report, other measures of the vaporization rate will be Introduced. 

Equation (4) gives the rate under normal vaporization conditions; It Is 
generally considered to be the maximum rate of vaporization. Steps or process- 
es that occur In series with the vaporization process can slow down the rate; 
but they can never speed It up. Only a process that parallels (bypasses) the 
normal vaporization step (e.g., chemical vaporization) can result In faster 
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vaporization. As we shall see, any of these steps may become the sole rate- 
controlling step. After we have completed the treatment of normal vaporiza- 
tion, we will discuss the other possible factors as rate-controlling steps: 
partial enclosure, matrix gas flow, chemical reactions, gaseous diffusion, 
presence of a coating, solid state diffusion - controlled Incongruent vaporiza- 
tion, and miscellaneous factors. 


Normal Vaporization 

In this section we will consider vaporization Into a vacuum (eq. (4)). At 
this time there Is no grounds for the selection of a back pressure, P(Y) (eq. 
(3)). However, In Part II of this report, the quantity P(Y) will be consid- 
ered because It arises out of the Interactions which are discussed there. 

Congruent vaporization . - When equation (4) Is used, the value obtained 
Is the rate of vaporization of a material from a clean surface under vacuum 
conditions. To use this equation, only the vapor pressure, and molecular 
weight have to be known. The answer obtained Is Independent of all parameters 
associated with the experimental setup. For a material such as Re which vapor- 
izes congruently, the calculated rate of vaporization will be virtually cons- 
tant as long as the solid material remains. 

The vapor pressure of rhenium Is 5.92xl0 -8 atm at 2773 K (ref. 3, p. 

447); Its molecular weight Is 186.2 g/mole. These values yield a normal vapor- 
ization rate, W(nv), at 2773 K of 6.8xl0 - ^ g/cm^ sec. In figure 2 the 
values of W are given as a function of temperature. Because the density of 
Re Is 20.53 g/cm J (ref. 4), the recession rate, G, Is 3.21xl0 -8 cm/sec or 
1.19 ym/hr. 

Incongruent vaporization . - Most structural materials are not pure sub- 
stances. Many are one phase alloys (solid state solutions). From a vaporiza- 
tion viewpoint, such one phase alloys behave In a similar manner to other 
solutions such as aqueous solutions and polymers containing plasticizers. In 
general when an alloy vaporizes, the vapor consists of more than one type of 
species. It Is likely that a different species Is present for each component 
of the alloy. Ideally for a single phase alloy, the partial vapor pressure, 
P(J)* for each component (j) would equal the vapor pressure of the component 
In Its pure form, P ' ( j ) , times Its mole fraction In the alloy, n(j). 

P(j) = n(J) P'(J) (7) 

The corresponding value of W(j) could then be found via equation (4) and the 
total rate of vaporization of all components would be the sum of the Individual 
rates. 

W = I W(J) (8) 

If all the components had the same vapor pressure and molecular weight, then 
all of the components would vaporize In the same ratio as they are present In 
the solid material. Such a condition would still be called congruent vaporiza- 
tion. However, In the usual case of an alloy, the pressures and the molecular 
weights of the components are not the same, so that the composition of the 
material would change as vaporization proceeds. This change In composition 
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(mole fractions of the components, n(j)'s) means that for Incongruent vapor- 
ization, the actual rate of vaporization changes with time and that the rate 
calculated via equations (4), (7), and ( 8 ) are only Initial rates when the 
Initial values of n(j) are used. 

Furthermore, most alloys are not Ideal solutions so that the values of 
W(j) obtained assuming Ideality are only approximations. If an experimental 
value for P(j) Is available, It should be used In preference to one obtained 
from equation (7). Note that In some alloys, only a single component may be 
volatile. This Is analogous to an aqueous salt solution where the water can 
vaporize leaving the salt behind. 

Some materials such as TaC also vaporize Incongruently even though they 
are generally considered to be compounds. No vapor species such as TaC Is 
known to exist: the vapor species over solid TaC are Ta, C, and Its polymers 

(e.g., the trlmer, C 3 ). Moreover, equation (7) Is not obeyed for either Ta 
or the C species. To understand better the vaporization of TaC, let us look 
at Its temperature - composition phase diagram In figure 3 (ref. 5, p. 84). 

At the one-to-one atom ratio, the TaC phase Is In equilibrium with solid carbon 
at 2773 K. This means that the partial vapor pressure of a carbon species 

arising from TaC Is the same as that for the respective one from solid carbon. 

Also, at the stoichiometric composition of TaC (ref. 5, p. 99), the carbon 
vaporizes much faster than the Ta, and the dominant vapor species Is the carbon 
trlmer, C 3 , because It has the highest vapor pressure above free carbon 

(ref. 3, p. 445). The second most Important species Is the carbon monomer, C. 

For C 3 , the value of the vapor pressure Is 1.48xl0 -4 atm at 2773 K. 
Equation (4) requires using the molecular weight of the trlmer (36.03 g/mole) 
for the value of M. Thus, the vaporization rate, W, Is calculated to be 
7. 48x1 0 -4 g/cm 2 sec. The vapor pressure of the C monomer species (M = 

12.01 g/mole) Is 7.17x10"^ atm, resulting In a W of only 2.09x10“^ g/cm 2 
sec. The value of W(nv) should be taken as the sum of W(C) and W(C 3 ), 
7.69xl0 - ^ g/cm 2 sec. However, In some portions of this report, the contri- 
bution of the C monomer Is disregarded because of the complications that 
would arise. In figure 2 are plotted the values of W for C and C 3 as 
a function of temperature. The graph shows that the vaporization rate of the 
C monomer species accounts for less than three percent of the total vaporiza- 
tion rate of the carbon from TaC. Here, we reemphasize that the calculated 
values for TaC are not expected to be constant because It vaporizes Incongru- 
ently. However, as long as the distinct carbon phase exists In the surface of 
the TaC material, the vaporization rate will equal the Initial rate. Subse- 
quent rates of vaporization are controlled by solid state diffusion which will 
be discussed much later. 

The vapor pressure of the Ta species over TaC Is about 9x10“^ atm 
(ref. 5, p. 101). From equations (4) and ( 6 ), one finds W(Ta) to be 
1.02xl0 -9 g/cm 2 sec and Z(Ta) to be 1.3xl0 -6 . This relatively low 
rate confirms the fact that the carbon vapor species are the dominant ones. 

When a material vaporizes Incongruently, the evaluation of G Is not 
always a straight forward calculation. Assumptions are sometimes required. 

If TaC Is a two phase region (as the phase diagram Indicates; fig. 3), then 
G can be obtained by using the density of graphite (2.25 g/cm^; ref. 4) In 
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equation (5). The resultant value of G Is 3.42x10“* cm/sec or 1.23x10* 
ym/hr. If one assumes that TaC consists of only a single phase and that Its 
density Is not very sensitive to carbon concentration, G Is obtained by using 
the density of TaC (13.9 g/cm 3 , ref. 4) In equation (5). The value of G 
from this calculation Is 5.53x10“ 5 cm/sec or 1.99xl0 3 ym/hr. A third and 
least likely assumption Involves the conversion of the TaC phase Into another 
phase (e.g., Ta) as the carbon vaporizes. The following equation would then 
be applicable. 


G = W (M(TaC)/D(TaC) - M(Ta)/D(Ta) ) / M(C) (9) 

The density of Ta Is 16.6 g/cm 3 . This third assumption gives a value for 
G of 1.89x10“* cm/sec which Is 6.82xl0 3 ym/hr. The phase relationships 
for the system under study determine which of these three methods should be 
used. Thus, for other types of Incongruently vaporizing materials, the last 
method of calculation may be the most preferable. However, when one is only 
Interested In order of magnitude values, no great care Is needed In the selec- 
tion of the equation. Throughout the rest of this report when G Is calcu- 
lated for TaC, we will use equation (5) with D equal to the density of 
graphite. 

In discussions Involving Incongruent vaporization, It Is helpful to use 
the term F(r), the fraction of the volatile component remaining In the solid 
material. This can be calculated as follows If W can be considered to be 
constant. 


F(r) = W T'/D E 1 S (10) 

In this equation, T', Is the time, S Is the original thickness of the materi- 
al, and E 1 Is the original weight fraction of vaporizing species. For the 
Incongruently vaporizing TaC, E' Is M(C)/M(TaC), 0.062. Equation (10) may 
also be used for congruently vaporizing material (e.g., Re); In this case, E' 
would be unity. Another concept that Is sometimes helpful Is the useful life- 
time, T". We have arbitrarily defined T" as the time at which F(r) = 0.5. 

If W is again assumed to be constant. 

T" = D E' S/2 W = D E' S/2 Z W(nv) (11) 

Equation (11) like equation (10) Is applicable to both congruently and Incon- 
gruently vaporizing material; this Is due to the presence of E' In these 
equations. The value of T"(nv) for TaC would be 5.6 sec while T"(nv) for 
Re would be 42 hr. Care must be taken when comparing lifetimes between vari- 
ous materials; many materials are unusable long before the value of F(r) 
would become reduced to 0.5. 

As far as vaporization Is concerned, our two model materials may be called 
extremes. While the vaporization of rhenium Is very small at 2773 K the vapor- 
ization of carbon component of the TaC Is relatively large. As we shall see 
later, chemical vaporization can even Increase the rate of vaporization for 
TaC. However, other factors may tend to decrease the actual rate of vaporiza- 
tion of TaC to a value less than the values calculated from equation (4). 

Presence of several phases . - In many alloys, the existence of multiple 
phases Is a necessity. The additional phases add strength or other desirable 
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qualities. For such a system, one must consider vaporization of each component 
(j) from each phase (1). One uses P(j,1) and M(J) via equation (4) to 
obtain the corresponding value of W(j,1). The value pf W(j), the overall 
vaporization rate of the jth component, should be obtainable as follows. 

W(j) = E(W(J,1) A(1)) / 2A(1) (12) 

where A(1) Is the area of the phase on the surface of the material. 

Equation (8) Is then used to calculate the total value of W. When the phases 
are In equilibrium, the calculations simplify because the vapor pressure of a 
specified component Is the same In all the phases. Equation (12) becomes un- 
necessary, and the calculations are the same as for a one phase system. How- 
ever, while the vapor pressure of the component Is the same In all 

phases. Its concentration Is not likely to be the same. Thus, as vaporization 
proceeds, displacement from equilibrium may occur. 

One can recall that at the stoichiometric composition of TaC (fig. 3), 
two phases exist: C and substolchlometrlc TaC. From the discussion In the 

previous paragraph, one can now see why the system was treatable as one phase. 
Another point to note from figure 3 Is that If equilibrium Is allowed to per- 
sist as vaporization occurs, the carbon phases will disappear, resulting In a 
truly one-phase system. 


Evacuated Partial Enclosure 

When vaporization occurs In an evacuated enclosure, the pressure closely 
approaches the vapor pressure of the material present; the rate of vaporization 
Is thus controlled by the area of the exit hole. The rate of vaporization per 
unit area of the exit hole Is given by equation (4). The average rate per unit 
Inside area of the enclosure Is obtainable by multiplying this value by the 
ratio of the area of the hole, A(h), to the Inside surface area, A(s). Thus 

W = a P(v) sqr(H) A(h)/A(s) = a P(v) R'/2 L (13) 

The final form of the equation Is applicable to our reslstojet where irR ' 2 

Is the hole area and irR'L Is the Inside area. It becomes obvious that 

Z = A(h)/A( s) = R'/2L = 0.0025 (14) 

and that Z Is Independent of the material used. Via equation (6) the respec- 

tive values of W for Re and TaC would be 1 . 7x1 0~ 9 and 1.87x10-° g/cm 2 
sec. Such average values will Increase with the size of the exit opening. In 
the present approximation, the vaporization Is assumed to be uniform over the 
entire Inside area. In a later section, we will see that when the existence of 
an enclosure Is coupled with the vaporization process, the values of W and Z 
are dependent on Y, the distance along the tube (enclosure). Equations (13) 
and (14) may be of use In many systems. However, for the reslstojet, they are 
applicable only during times that the exit tube may be at temperature without 
the presence of a matrix gas. 

One should note that vaporization from an evacuated container with an 
orifice (small properly shaped hole) Is the basis of the Knudsen method of 
determining vapor pressures. The vapor pressures are calculated, of course, 
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via equation (4). Also note that In Part II of this paper, discussion of the 
Interaction of this mode of vaporization will be minimal because this mode can 
Interact with very few of the other modes. The other modes, In general, re- 
quire the presence of a matrix gas while this mode requires It absence. 


Flowing Matrix Gas 

The first step In the treatment of the effect of a flowing matrix gas In 
an enclosure Is to express Its flux, F(m), In terms of Its flow, F'. 

F(m) = b T' M(m) P(m)/A(h) (15) 

where m refers to the matrix gas, and b Is 1/RT and numerically equal to 
4.4x10“° cm- 3 atm- 1 (4.4xl0“ 12 ergs" 1 ) at 2773 K. Next, the flux of 
the vaporized species, F(Y), can be written In terms of the flux of the matrix 
gas. 


F(Y) = F(m) P(Y) M/P(m) M(m) = b F 1 H P(Y)/A(h) (16) 

where Y Indicates the vaporized species. The expression for W can be ob- 
tained by setting It equal to the flux of the vaporized species parallel to 
the flow multiplied by the ratio of the exit hole area to Inside surface area. 


W = F ( Y ) A(h)/A(s) 


= b F' P(v) M/A( s) = 3.75X10' 4 P(v) M 
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The second form equation (17) results from combining Its first form with equa- 
tion (16). The final form of equation (17) holds for our reslstojet design. 
Combining equations (4), (6), and (17) yields 


Z = b F‘ sqr(M)/a A(s) = 1.66xl0“ 6 F 1 sqr(M) = 4.45x10 4 sqr(M) (18) 

The W and Z given by equations (17) and (18) are average values. When 
we discuss the coupling of matrix gas flow with normal vaporization, we shall 
see that W and Z are dependent upon the distance along the tube or enclo- 
sure. From equation (17) the values of W for Re and TaC were found to be, 

respectively, 4.13xl0 -9 and 2.0xl0“& g/cm 2 sec. From equation (18), the 
respective values of Z are 0.00607 and 0.00267. Thus, the effect of the ma- 
trix gas at the prescribed velocity (for the reslstojet) on the average rate of 
vaporization of the tube material Is very roughly the same as If there were no 

matrix gas present: the values of Z are close to 0.0025 (eq. (14)). From 

equation (18), one can see that decreasing the flow rate, F 1 , will decrease 
Z, and visa versa. Thus, with a slow enough flow, the value of Z could be- 
come quite small. 


In the extreme where F Is zero, equation (18) states that Z Is also 
zero. This Incorrect answer from equation (18) Is the result of approximations 
used Its derivation. However, closer approximations need not be used because 
the situation Is trivial In that the matrix gas would be quickly lost Into the 
surrounding vacuum. After the matrix gas Is lost, equation (14) gives the 
proper steady state value for Z. Another way to look at this problem Is given 
In the following section. 
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Stationary Matrix Gas 


To have a stationary (nonflowing) matrix gas, the external pressure of 
the matrix gas (that outside the tube) must be equal to Its pressure within 
the tube. Of course, this condition Is not possible for the reslstojet, but 
It may exist for other design problems. The logical way to attack such a 
problem Is via gaseous diffusion. However, until we discuss the coupling of 
steps In Part II, the equation for this type of diffusion process cannot be 
set up. For now, we will be satisfied to calculate the maximum values of W 
and Z for this process. These are the same values given for an evacuated 
enclosure (eqs. (13) and (14)). Thus for all materials, Z(max) = 0.0025 when 
a stationary matrix gas Is present; W(Re) = 1 . 7x1 0~ 9 g/cm 2 sec, and 
W(TaC) = 1.87X10- 6 g/cm 2 sec. 


Chemical Vaporization 

The presence of a matrix gas can also give rise to the possibility of a 
chemical reaction between the material and the gas. Itself. This could result 
In an Increased rate of vaporization because a volatile species may be formed 
during the chemical reaction. For the metal, rhenium, a chemical reaction Is 
very remote; Re should be unreactive chemically under the usual conditions of 
the reslstojet. The* situation with tantalum carbide Is quite different. For 
the time being, the carrier (matrix) gas Is assumed to be hydrogen, and the 
reaction product Is assumed to be acetylene (C 2 H 2 ). This Is the most like- 
ly product from equilibrium considerations (ref. 6). Because of the excess of 
carbon In stoichiometric TaC, It Is valid to write the chemical reaction as 
follows. 


2C(sol1d) + H 2 (gas)^ZZ^C 2 H 2 (gas) (19) 

The opposing arrows Indicate that both forward and reverse reactions can occur. 
The rate of the forward reaction should be proportional to the pressure that 
the hydrogen exerts (called the hydrogen partial pressure), and rate of the 
reverse reaction Is proportional to the partial pressure of the C 2 H 2 . 

When reaction (19) Is at equilibrium, the rate of the reverse reaction becomes 
equal to the rate of the forward reaction so that one may write 

K(p) = P(C 2 H 2 )/P(H 2 ) (20) 

where K(p) Is called the equilibrium constant for the reaction. Many such 
constants are tabulated In reference 6 and given In diverse sources In the 
chemical literature. At 2773 K, the log K(p) (base 10) for reaction (19) Is 
equal -1.453 resulting In a partial pressure for C 2 H 2 of 0.1113 atm when 
the partial pressure of hydrogen pressure Is 3.16 atm. When chemical reac- 
tions are Involved, equations (3) and (5) must be modified respectively as 
follows. 


W(cs) = a (P(cs) - P( Y,cs) ) M'/sqr(M) (21) 

G(cs) = a (P(cs) - P(Y,cs)) M'/D sqr(M) (22) 
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where P(cs) Is the vapor pressure of the chemical species, and P(Y,cs) Is 
the pressure of this chemical species In contact with the material. In the 
present discussion, P(Y,cs) = 0. The M' Is the molecular weight of the sub- 
stance lost by vaporization, while M Is the molecular weight of the actual 
vapor species which Includes a weight contribution from the matrix gas (e.g., 
h 2 ) due to the chemical reaction. When considering the formation of 
C 2 H 2 , H Is the molecular weight of C 2 H 2 (26.036 g/mole), and M' Is 
the molecular weight of C 2 (24.02 g/mole). Thus, W Is 0.441 g/cm 2 sec, 
the Initial rate of vaporization of carbon from TaC at 2773 K via chemical 
reaction with hydrogen to form acetylene. As defined In equation ( 6 ), the 
value of Z Is equal to the vaporization rate of the chemically formed spe- 
cies (cs) divided by the rate for the normal vaporizing species (nv). Thus, 

Z = P(cs) M'(cs)/P(nv) sqr(M(cs) M(nv)) (23) 

In the present case, "nv" refers essentially to C 3 while "cs" refers to 
C 2 H 2 with H'(cs) = M(C 2 ). Examination of equation (23) Indicates that 
when a chemical reaction Is Involved, It Is simpler to calculate Z(cs) by ■ 
dividing the numerical value of W(cs) by the numerical value of W(nv), 
7.69x10-4 g/cm 2 sec. Thus, It follows that Z = 574; this shows that the 
vaporization rate by chemical reaction could be much greater than that by the 
normal physical process. Via equation (22) and the assumption that the loss 
In weight Is from the carbon phase, the value of G Is calculated to be 0.197 
cm/sec. Remember that W, G, and Z are maximum values. If equilibrium Is 
not reached, the actual rate could be much less, even less than that for sim- 
ple vaporization. The only way to determine the real kinetic value Is via 
experiment. The value of T" from equation (11) Is 9.7xl0 -3 sec. 

The TaC may also react with H 2 to form methane, CH 4 . 

C(solld) + 2H 2 (gas) V-^-CH 4 (gas) (24) 

At 2773 K, the log KtpUtnr 1 )) for reaction (24) Is -4.076 (ref. 6 ). For 
this case, the rate of the forward reaction should be proportional to the 
square of the hydrogen pressure because two hydrogen molecules are Involved In 
the reaction. Thus, 


K(P) = P(CH 4 )/P(H 2 ) 2 (25) 

The value of P(CH 4 ) Is therefore equal to 8.38xl0 -4 atm. Using equa- 
tions (5) and (21) with M' = M(C) and M = M(CH 4 ), one obtains U = 

0.00212 g/cm 2 sec and Z = 2.8. Thus, at the stoichiometric composition of 
TaC, the vaporization via methane formation Is expected to be slower than that 
by acetylene formation but greater than that for C 3 . The value of Z(cs) 
would be 577, the sum of the Z for C 2 H 2 and the Z for CH 4 . 

When TaC Is exposed to a partial pressure of oxygen, another chemical 

reaction can occur. The reaction can be written 

C + 0 2 -^Z?r C0 2 (26) 
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The value of K(p) for reaction (26) Is 2.92xl0 -7 . The value of P(C02) 
could be calculated In the same manner as we did for acetylene. However, this 
calculation can be bypassed because K(p) Is very large Indicating that the 
equilibrium Is far to the right; this means that any oxygen present Is almost 
completely converted to CO2 as long as enough carbon species Is present. 

At equilibrium, the partial pressure of CO2 would be equal to the partial 

pressure of the oxygen originally In the system. The value of W can then be 

calculated from equation (21) by using this assumption.. For an assumed Ini- 
tial oxygen partial pressure of 0.1 atm, the value of W Is 0.152 g/cm* 
sec. The corresponding value of Z Is 198. Note that In this case, the 
exact value of K(p) need not be known. 

The presence of oxygen In a matrix gas of hydrogen gives rise to the fol- 
lowing reaction. 

0 2 + 2H 2 2H 2 0 (27) 

The equilibrium for this reaction Is far to the right; K(p) Is 50.75 recipro- 
cal atm. Thus, before the oxygen contacts the TaC, It can react with the hy- 

drogen matrix gas, and the resultant water vapor can react with the TaC to 
form CO2 along with the recovery of the hydrogen which originally reacted 
with the oxygen. This can give the same result as If the oxygen reacted di- 
rectly with the TaC. The final partial pressure of the H2O Is many times 
greater than that of the O2, but both of these are negligible when compared 
with the partial pressure of the CO2. The existence of reaction (27), 
therefore, does not nullify any calculations made for reaction (26). Moreover, 
If one starts with H2O vapor, the H2O can also react with the TaC to form 
CO2. Because there Is only one oxygen atom In H2O compared with two In 
CO2, an original partial pressure of 0.1 atm H2O would result In a final 
0.05 atm of CO2. This would result In half the rate of vaporization that 
would have resulted from a 0.1 atm of oxygen. 

The other refractory carbides such as those of hafnium, titanium, and 
tungsten would be expected to react In a similar manner with hydrogen or water 
vapor. In fact for those of hafnium and titanium, the activity of carbon at 
the stoichiometric composition Is unity just as It Is In TaC. Therefore, the 
value of W(cs) for HfC (eq. (21)) would be the same as to TaC. Chemical 
vaporization Is also a common phenomenon for the refractory metals. Metals 
such as tungsten, molybdenum, and even platinum will react with oxygen and 
vaporize chemically In the form of their respective volatile oxides. Experi- 
mentally, the rate of oxidative vaporization has been found to vary with the 
oxygen pressure as one would predict from the equilibrium considerations. 

These are only a few of the cases where materials vaporize chemically. To 
decide If a particular material Is going to vaporize chemically, one should 
first ascertain what gases are present and then from equilibrium data deter- 
mine to what degree reaction will occur. Of course In high vacuum or with 
only noble gases present, chemical vaporization Is minimized. 

Note that vaporization (when It Is controlled exclusively by a chemical 
reaction step) Is again not affected by parameters of the experimental setup, 
such as the length of the tube. Its maximum rate Is determined by the 
thermodynamic data for the reaction that Is occurring. The actual value Is 
determined by the rate of the reaction; this would have to be found 
experimentally. 
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In addition to leading to the formation of a gas species, a chemical re- 
action can cause the formation of a solid phase. This Is a very common occur- 
rence for many metals such as Iron. Typically, such a surface phase Increases 
In thickness with time. If the surface phase Is volatile, the calculation of 
Its vaporization rate would be given by equation (4) with the vapor pressure 
used being that of the surface phase material. Vaporization of the substrate 
(the phase underneath) can also occur via solid state diffusion through the 
surface phase. This Is discussed in a later section called "Surface Coating". 


Diffusion Through a Gaseous Boundary Layer 

Another way for a matrix gas to affect vaporization Is by what we shall 
call back reflection. Vaporizing molecules will collide with molecules of the 
matrix gas, resulting In a change in their direction. Some of these vapor- 
izing molecules will thus be directed back to the solid material and condense. 
One can also look at the problem in terms of a vaporizing species diffusing 
through the matrix gas. It Is through diffusion theory that the problem can 
most easily be solved. 


In a relatively rapid flow of the gas, one may assume the presence of a 
boundary layer through which the vaporizing species must move In order to leave 
the surface of the tube material and get to the main stream of the flowing 
matrix gas. The actual flow may be either laminar or turbulent, but at pres- 
ent we need not concern ourselves as to which It Is; all that Is necessary Is 
the assumption of the existence of a boundary layer. With our symbolism, the 
diffusion equation (ref. 7) is written as follows: 


W = - 0(g) dC"(v)/dX = - b M D(g) dP/dX 
= b D(g) M P/S' 


(28) 


C"(v) Is the carbon concentration In the gas phase In g/cm 3 ; X Is the dis- 
tance out from the surface of the solid material; W Is the flux In that direc- 
tion; D(g) Is the gaseous diffusion coefficient In cm 2 /sec; and dC"(v)/dX 
Is the concentration gradient of the diffusing gas species perpendicular to 
the surface. The second form of the equation was obtained by using the Ideal 
gas equation with b = 1/RT. The final form of the equation assumes planar 
diffusion. This assumption can be made because the thickness of the boundary 
layer. S' (cm), should be small compared with the radius of the tube, R'. Also, 
In the final form of equation (28), we assume that the partial pressure (con- 
centration) of the vaporized (diffusing) species Is zero In the bulk of the 
matrix gas, and that steady state conditions are achieved. One may use the 
vapor pressure for the value of P because this Is the maximum value that P 
can have. The equation for Z Is obtained by combining equations (4), (6), 
and (28). 

Z = b D(g) sqr(M)/a S' = 5.21xl0" 6 D(g) sqr(M)/S' (29) 

The value of D(g) required In equations (28) and (29) can be calculated 
from the following equation (ref. 2, pp. 268 and 281). 


14 


D(g) = sqr(V(v) 2 + V(m) 2 )/3 a N" Q ,,2 ( 1 + A) 

= (2RT/ir) 3/2 sqr(l/H(v) + l/M(m))/3 N P Q" 2 

= 3.07X10 -14 sqr(l/M(v) + 1/M(m))/P Q" 2 (30) 

where the V's are average molecular speeds; v refers to the vaporizing spe- 
cies and m, to the matrix gas; N" Is the number of molecules per cm 3 ; and 
Q" Is the effective average molecular diameter (cm). The A' Is a term to 
correct for persistence of direction after collision and Is dependent on the 
Interaction between molecules; having no simple way to evaluate A', we Ignore 
It. In the final forms of the equation, we have made proper substitution for 
V's, etc. The N Is Avogadro's number, and P Is the total pressure In atm. 
If we select a reasonable value for Q" (4xl0” 8 cm, ref. 2, p. 285), the 
values of D(g) for Re and TaC are, respectively, 4.3 and 4.4 cm 2 /sec at 
2773 K and 3.16 atm. It Is also possible to estimate the value of D(g) from 
the known experimental data of similar molecules (e.g., ethane; Its D(g) val- 
ue Is 5.12 cm 2 /sec at 2773 K and 3.16 atm of hydrogen gas; ref. 8). Thus, 
our calculated value for 8(g), 4.4 cm 2 /sec, Is deemed to serve as an approx- 
imate number that can be used later for the value of D(g) for other similar 
carbon-containing gaseous species. 

Using the calculated values of D(g), we obtained the curves In figure 4, 
plots of W as a function of the boundary layer thickness, S'. From this 
figure, one sees that W and Z decrease with Increasing boundary layer 
thickness. If one would assume too small of a value for S', the calculated 
value of Z would be greater than unity. This would represent an unaccept- 
able physical situation. A reasonable assumption for the boundary layer Is 10 
percent of the diameter; this means that S' Is 0.005 cm. The respective 
values of W for Re and TaC would be 4.16xl0” 8 and 2.07xl0~ 5 g/cm 2 sec. 

The respective values for Z are 0.061 and 0.028. The existence of a bound- 
ary layer In the reslstojet thus appears to reduce only moderately the vapor- 
ization rates for Re and TaC. Again, calculations are Independent of the 
parameters of the experimental setup. The Important parameters are the pres- 
sure of the matrix gas, the thickness of the boundary layer, and the gas dif- 
fusion coefficient. 


Surface Coating 

In many practical systems, materials are coated with some protective layer 
to Inhibit heat transfer or to shield the substrate from oxidation. However, 
no such coating Is planned for the reslstojet. As discussed before, coatings 
may also exist due to formation via chemical reaction between the solid and 
the matrix gas. 

The method of calculating the vaporization rate from a coating layer would 
again be given by equation (4) by using the proper values of P(v) and M. 
However here, we are more Interested In the vaporization of the material under- 
neath the coating (the substrate). A substrate may vaporize even though cov- 
ered with a protective coating; It Is possible for the volatile species In the 
substrate to diffuse through the coating and then evaporate. The solution of 
this problem requires the use of the steady state diffusion equation. The 
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steady state equation for solid state diffusion has the same differential form 
as that for gaseous diffusion (eq. (28)). 

W = - D(s) dC"(s)/dX = D(s) E' D/S" (31) 

where dC"(s)/dX Is the gradient of the diffusing (vaporizing) species and 
D(s) Is Its diffusion coefficient In the coating. The final form of the 
equation assumes planar diffusion and the gradient to be uniform across the 
coating with the concentration of the diffusing species equal to zero at the 
coating-gaseous Interface. D Is the density of the coating (g/cm 3 ), E' Is 
the solubility (g/g) of the dlffuslng-vaporlzlng species In the coating, and 
S" Is the thickness of the coating. The value of Z Is obtained by dividing 
equation (31) by equation (6). 

Z = D(s) E' D/a P(v) S" sqr(M) (32) 

Any evaluation of equations (31) or (32) for our reslstojet design will 
be extremely approximate because at the present time no coating Is contem- 
plated for the system. To apply equations (31) and (32) to the reslstojet, we 
will assume 10 ppm solubility (E' = 0.00001), a density of 10 g/cm 3 for the 
coating, and a D(s) value of l.llxlO" 7 cm 2 /sec for both Re and TaC. A 
low solubility Is required for a coating to Inhibit vaporization. This value 
of D(s) correspond* to the value for carbon diffusion In TaC (ref. 9). By 
using these values and equation (31), U Is plotted as a function of S" In 
figure 5. For a reasonable thickness (0.01 cm), the value of U Is about 
l.lxlO -9 g/cm 2 sec for both Re and TaC. The fact that W for TaC and Re 
are the same In our calculation Illustrates that diffusion coating-controlled 
vaporization can be totally dependent on the diffusion parameters and Indepen- 
dent of the vapor pressure of the substrate. Remember, we have assumed the 
values of the diffusion parameters; their actual values are, of course, ex- 
pected to depend on the nature of the substrate. From equation (32) and our 
assumptions, the values of Z are found to be 1.6xl0 -3 for Re and 
1. 43X10- 6 for TaC. 


Solid State Diffusion-Controlled Incongruent Vaporization 

Because the title to this section Is quite long, we will refer to It In 
the future as dlffuslon/lncongruent vaporization. 

For rhenium, as well as other elements, the question of Incongruent vapor- 
ization does not arise because It contains only one element. Even many com- 
pounds vaporized congruently. But as we have seen before, TaC vaporizes 
Incongruently: the elements present do not vaporize at the same rate. The 

same Is true for the majority of alloys and structural materials. In TaC, the 
carbon vaporizes first and Ta Is left behind. As this happens In a one phase 
system, the partial vapor pressure of C decreases and so does Its vaporiza- 
tion rate. The same sort of thing happens In other Incongruently vaporizing 
systems. In these cases, vaporization can become controlled by the solid 
state diffusion of the volatile species to the surface. Thus, the time 
dependent form of the diffusion differential equation becomes applicable (ref. 
7, p. 9). 
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(33) 


dC"(s)/dT' = D(s) d 2 C"(s)/dX 2 

where C"(s) Is the concentration of the diffusing species. 

The solution of this differential equation for actual systems Is complex 
because the vapor pressure of the vaporizing species (carbon) Is not necessar- 
ily a well behaved function of Its concentration. Besides, all of the experi- 
mental data are not available. It Is therefore necessary to make some 
assumptions, some of which may not be completely true. For example In the TaC 
system, we will assume that a single phase exists between the compositions TaC 
and Ta. If solid state diffusion Is to be the sole rate-controlling step, one 
must assume that the volatile component vaporizes as soon as It reaches the 
surface layer. This Is the same as assuming an essentially Infinite vapor 
pressure and the surface concentration of the volatile species to be zero. 

With this as a boundary condition, the solid state diffusion equation (eq. 
(33)) may be Integrated (ref. 7, p. 45). By using our symbols and assuming 
planar diffusion, 


F(r) = 8 E (exp(-D(s)(2n + l) 2 ir 2 T'/4S 2 )/(2n + l) 2 ir 2 ) (34) 

n=0 

where F(r) Is the average bulk carbon concentration, F(r) being defined as 
unity for the stoichiometric TaC composition. 


W = 20 D(s) M(C) £ (exp(-D(s)(2n + l) 2 ir 2 T'/4S 2 )/M(TaC) S) (35) 

n=0 


The computer program for equations (34) and (35) are given In table 1(a). The 
value of Z can be calculated via equation (6). The present model Is re- 
ferred to as Model 0 to distinguish It from models presented later which In- 
clude other vaporization steps. 

In figure 6, a plot of F(r) versus the square root of time shows para- 
bolic dependence. At F(r) = 0.5, T 1 Is 168 sec so that 168 sec may be taken 
as the value of the half life, T". In figure 7, the relative vaporization 
rate Is given as a function of time. The Initial value of Z Is, of course, 
unity, meaning that W = W(nv). The final value of Z Is obviously zero. 

This Is a result of the depletion of carbon In the TaC, leaving Ta behind. 
Because Z Is a function of time, T', for Incongruent vaporization. It Is 
convenient to define, Z', pseudo Z for an Incongruent process which would be 

Independent of time. The definition clearly requires a comparison of some 

time Independent rate parameter with the corresponding one under the standard 
vaporization conditions. We chose to define V as follows. 

V = T"(nv)/T" (36) 

The value of V for the Incongruent vaporization of TaC would equal 0.033. 

Many times parabolic kinetics as noted In figure 6 means the formation of 

a new phase on the surface of the material as In the oxidation of some metals. 
Even though a postulate of Model 0 Is the existence of a single phase, we 
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thought It would be wise to show that Model 0 does not predict a sharp demar- 
cation between a carbon-rich and a carbon-poor region. This is done in figure 
8 where the carbon concentration profiles are plotted for various values of 
F(r). Again, an Infinite series equation was required (ref. 7, p. 44). In 
figure 8, there Is Indeed no suggestion of such a sharp concentration 
demarcation. 

According to equation (34), F(r) can be expressed as a function of 
T'/S 2 . This means that T'/S 2 Is a constant for any given value of 
F(r). It follows from this that the half life, T", Is proportional to the 
square of the tube thickness. Thus, a slight Increase In the thickness of the 
reslstojet tube can add significantly to the value of T". 

Model 0 Is, Indeed, an oversimplified solution to the problem. But It 
will probably give a correct order of magnitude answer for TaC. The total 
solution would Involve the diffusion of C through Ta 2 C as well as through 
TaC (fig. 3). Treatment of the problem to this mathematical depth would re- 
quire a knowledge of the exact dependence of carbon activity on Its concentra- 
tion. Furthermore from a practical viewpoint, when the composition of Ta 2 C 
Is approached, the life of the material Is essentially at an end, and calcula- 
tion past this point would be only of academic Interest. For other materials, 
a similar logic should hold so that a one phase diffusion model Is sufficient 
for many practical problems. 

Model 0 Is applicable to many space problems where an alloy at an ele- 
vated temperature Is exposed to the vacuum of space. The only unique parameter 
needed for the application of Model 0 Is the diffusion coefficient of the 
vaporizing species In the alloy. If It Is unknown. It can be approximated. 

Note that for Model 0, the partial vapor pressure of the vaporizing species 
does not have to be known. The F(r) calculated Is the average fraction of 
the vaporizing species remaining In the alloy. Model 0 Is the most severe 
model one can use. Models 1 and 2 to be described later requires knowledge of 
the vapor pressures and should result In a longer predicted life. 


Miscellaneous Factors 

There are still other factors that can affect the rate of vaporization of 
a material: adsorbed gas on the surface, roughness of the surface, grain size, 

and surface tension. Surface tension effects can of course be calculated, but 
these are second order effects on vaporization. Thus, they are not Important 
In the present discussion. In some cases, other miscellaneous effects may 
have been Investigated for a particular material. However, a general treatment 
of such effects Is not available, and the effect of such factors (variables) 

Is expected to be quite small relative to the effects discussed In this report. 

Degassing of adsorbed gases may have a notable effect on the physical 
properties of some materials. However, we desire to keep this effect outside 
the scope of this report. 


18 


PART II: COUPLING OF KINETIC STEPS 


General Principles 

Because vaporization Is a kinetic process, It couples with other kinetic 
steps In a predictable manner. When the process Is Y and T' dependent, 
the process must be treated as a special case. However, when the process Is 
Y and T' Independent, general equations can be used as one will see In this 
present section. 

For several kinetic steps that work In parallel (l.e., Independently of 
one another using separate paths), the total rate Is equal to the sum of the 
Individual rates, W(k). 

W( total) = l W(k) (37) 

It follows that 

Z(total) = l Z(k) (38) 

The total rate Is greater than any Individual step. Also, If the rate of one 
step Is very much faster than the others, the overall rate Is controlled by 
the fastest rate. 

A similar equation describes the effective rate, W(ef), for the combina- 
tion (coupling) of several series steps If their kinetics are first order. 

W(ef) = 1 /l (1/W(k) ) (39) 

Similarly, 

Z(ef) = 1/2 (1/Z(k)) (40) 

A first order step Is one for which the rate Is dependent on the first power 
of the concentration (activity) of the species Involved. All the vaporization 
steps for Re are first order, but the normal vaporization step for TaC Is ex- 
pected to be third order as will be discussed later. For the combination of 
series steps, the overall rate Is slower than the slowest Individual step In- 
volved. If one of the steps Is very much slower than the others, the overall 
rate Is essentially controlled by the slowest step (ss). 

W(ef ) = W(ss) (41) 

Z(ef ) = Z(ss) (42) 

Note that equations (39) and (40) are symmetrical with respect to the various 
steps which are occurring. Thus, to use equations (39) and (40), one does not 
need to know the chronology of the steps. 

Equation (39) Is easily shown to be true for a two step process. Let us 
call W(1 ) and W(2) the rates of two kinetic processes when they are occur- 
ring by themselves. When they occur In series and steady state Is attained, 
the rates of both steps must be equal. Also, the rate of the first step Is 
given as follows: 
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( 43 ) 


W(ef) = (1 - C') W(1) 
The rate of the second step Is 


W(ef) = C' W( 2) (44) 

where the 1 In the factor (1 - C') of equation (43) represents unit activity 
(concentration) at the start of the first step. This first step reduces the 
activity of the volatile component to C' on the surface of the material. 

The value C' will contribute to the reverse reaction so that the net forward 
reaction Is proportional to (1 - C'). The use of this factor Is similar to 
the use of the factor, P(v) - P(Y) (eq. (3)). The second step starts with the 
volatile component at an activity of C' and takes It to zero; thus, W(2) Is 
multiplied by C', the difference between C' and zero. Elimination of C' 
between equations (43) and (44) yields 

W(ef) = W(l) W(2)/(W(1 ) + W( 2 ) ) 

= 1/(1/W(1) + 1 /W( 2) ) (45) 

which Is equation (39) applied to a two-step reaction. Note that equation 
(45) holds even If the system Is not Ideal. For a nonideal system, C' would 

equal the activity but not the concentration of carbon. In spite of this, 

equations (43) and (44) could still be combined In the same fashion. 

When a step In the overall vaporization process Is not first order, equa- 
tion (45) may not be valid. The most likely step not to be first order Is the 

normal vaporization step, Itself. In the case of TaC, three carbon atom In 
the solid must combine to form the trlmer carbon vapor species. 

3C(sol1d) ^ - C3(vapor) (46) 

Thus, the vaporization rate should be proportional to the cube of the carbon 
activity on the surface In the solid; such kinetics are referred to as third 
order. If the second rate-controlling step occurs after the third-order step, 
the C 1 In both equations (43) and (44) becomes equal to the cube of the car- 
bon concentration, and equation (45) would still be valid. However, If the 
third-order step comes after the other step, equation (44) must be modified as 
follows. 


W(ef) = C' 3 W(C 3 ) + C' W(C) (47) 

We have Included In equation (47) the contribution from the vaporization of 
the carbon monomer. By eliminating C' between equations (43) and (47) one 
obtains 


1/W(ef ) = 1/W(1) + 1/(W(C) + W(C 3 ) (1 + W(ef)/W(l)) 2 ) (48) 

The equation requires use of successive approximation, but W Is, neverthe- 
less, readily calculated. The Z Is given by a similar equation. 

1/Z(ef) = 1/Z(1) + 1/(Z(C) + Z(C 3 ) (1 + Z(ef)/Z(l)) 2 ) (49) 
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In contrast with equations (39) and (40), equations (48) and (49) are not sym- 
metric with respect to the steps involved. Equations (48) and (49) apply to 
the vaporization of TaC when the first step occurs before the vaporization 
step and the second step Is either the vaporization step or one occurring af- 
ter It. However, when one of the steps Is much slower than the other, equa- 
tions (41) and (42) can be used In place of equations (48) and (49). 

To prevent our subsequent discussion from becoming, too verbose, we have: 
limited our discussion to the coupling of two steps; used equations (41) and 
(42), where applicable; confined calculations to values of Z, then possible; 
and given much attention to the order In which the calculations are presented. 


Coupling Normal Vaporization With Other Steps 

Couplings Independent of Y and T‘ . - From Part I of this report, one 
can note that Z « 1 for all steps except for normal vaporization and chemi- 
cal vaporization. Therefore when dependence on Y and T' Is Ignored, the 
rate (Z) of normal vaporization coupled with one of the slow steps (ss) Is 
given by equation (42). This equation applies to the coupling of normal va- 
porization with an evacuated enclosure, flowing gas, stationary gas, boundary 
layer, or surface coating. These values of Z have already calculated In 
Part I of this report. Remember that many of these values are dependent on 
the parameters of the design and should be recalculated when another design Is 
Involved. 

Before we leave the subject of Y and T' Independence, we thought It 
wise to derive the effective value of Z by a formal route at least once. We 
have chosen to do this for the evacuated enclosure. Starting with equation 
(39), one Inserts the analytical expressions for the W's from equations (4) 
and (13). The resultant equation follows: 

W(j,1) = a P(j,1) sqr(M( j ) ) A(h)/(A(s) + A(h)) (50) 

Equation (50) Is written In the form for a multicomponent, multiphases system. 
The total contribution of the j*h component from all the phases Is obtained 
by using equation (11) and the total contribution of all the components can be 
computed via equation (7). But for our model materials, Re and TaC, we need 
consider only one component and one phase, so that we can drop the 1 and j 
In equation (50). Note that In subsequent equations where 1 and J are not 
specified, one merely has to Insert the 1 and j Inside to parentheses In 
order to consider a multicomponent and multiphase system. 

Equation (50) without the j and 1 Is not much different than equation 
(13); for the reslstojet with a very small hole area with respect to the In- 
side area, the correction Introduced by equation (50) Is Insignificant. How- 
ever, when the hole area Is relatively significant, equation (50) should be 
used. 

It follows from equation (50) that W for Re Is 1.7xl0 - ^ g/cm^ sec 
and W for TaC Is 1.87xl0~ 6 g/cm‘ sec. From these numbers and equation 
(5), Z for Re Is found to be 0.0025 and Z for TaC Is also 0.0025. These 
are the same value calculated when the effect of the evacuated enclosure was 
assumed to be the sole rate-controlling step (eq. (14)). 
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For the cases of the evacuated enclosure, flowing gas, and stationary 
gas, the values of Z obtained by using equation (42) should be referred to 
as average values, Z(av), because the actual Z values can be functions of 
Y. From a practical point of view, we are Interested In only one Y posi- 
tion, the place where vaporization Is the greatest because that position 
determines the life of the entire tube. However, first one must determine the 
value of Y for which the vaporization rate Is a maximum. 

Dependence on Y for evacuated enclosure . - If the pressure caused by 
the vaporizing gas were high so that there would be considerable collisions 
between molecules, one could treat the present problem with diffusion theory. 
However, at the low pressures involved, a vaporized molecule will either exit 
through the hole or collide with the wall of the tube somewhere. Starting 
with this model and the concept that the direction of a molecule Is random, 
the fraction of vaporizing molecules that exit through the hole Is equal to 
the solid angle that the hole makes with a single molecule divided by the to- 
tal solid angle. Thus, 

W(Y) = W (1 - cos(arctan(R'/(L - Y))))/2 

= a Z P( v) sqr (M)(l - cos(arctan(R/(L - Y))))/2 (51) 

where W and Z refer to the steps which are being coupled with the enclo- 
sure step. It follows that 

Z(Y) = Z (1 - cos(arctan(R'/(L - Y))))/2 (52) 

When one Is dealing with normal vaporization, Z = 1. Equations (51) and (52) 
give the rate of vaporization per unit area of the Inside surface of the tube 
as a function of Y, the distance from the entrance. The value of Z(Y) Is 
plotted as a function of Y In figure 9. One can see that most of the vapor- 
ization occurs near the exit of the tube. Thus, the lifetime of the tube Is 
determined by the value of Z(Y) at the exit (0.5). The value of Z(Y) Is 
greatest at the exit because the molecules near the exit possess a larger sol- 
id angle of escape. At a position 1 cm from the exit the value has dropped 
over two orders from Its value at the exit; Z(entrance) = 6.2xl0~&. By sum- 
ming Z over the length L and dividing by L, we found the average value 
of Z to be 0.00249 as compared with 0.0025 given In equation (14). If the 
enclosure were not circular or uniform, equations (51) and (52) are still ap- 
plicable with R' being the effective radius at the exit. 

Flowing matrix gas . - The coupling of a vaporization step with the flow- 
ing matrix gas step can also be considered to be a function of Y. At the 
entrance of the tube, we can postulate that the matrix gas Is free of vapor- 
ized material. It follows that this Is also the location where the rate of 
vaporization Is expected to be the greatest. At the exit hole, the matrix gas 
is expected to be nearly saturated with the volatile species, and so here the 
rate of vaporization should be the least. No longer Is equation (17) valid: 

W Is no longer proportional to F(Y). Instead, the value of W should be 
added to the value of F(Y) as one proceeds from the entrance to the exit of 
the reslstojet tube In the following manner. 
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dF( Y) = W A(surface area of a dY thick s!1ce)/A(cross section) 
or 


dF( Y)/dY = 2 W ir R'/ir R' 2 = 2W/R' 

= 2 a Z (P(v) - P( Y) ) sqr(H)/R' (53) 

where W and Z In equation (53) refer to the step being coupled with the 
flowing gas step, dF(Y)/dY Is the gradient of the vaporized species flux 
parallel to the length of the tube, and Y Is the distance from the entrance 
of the tube. Note that P(Y) can no longer be considered to be zero but Is 
the partial pressure of the vaporized species In the matrix gas. Moreover, It 
Is a function of Y. Equation (16) Is again valid for F(Y). Combining equa- 
tion (16) and (53), followed by Integration, gives an expression for P(Y). 

Then via equation (3), one can obtain 

W(Y) = a Z P(v) sqr(M) exp(-Y R'Z sqr(2ir/bM)/F' ) 

= 0.842 Z P( v) sqr(M) exp(-225 Y Z/sqr(M)) (54) 

The corresponding value of Z(Y) would be the exponential factor times Z. 

Z(Y) = Z exp(-Y R'Z sqr(2ir/bM)/F ' ) 

= Z exp(-225 Y Z/sqr(M) ) (55) 

When the step being coupled to the flowing gas step Is normal vaporiza- 
tion, Z = 1. By using the values of the parameters given before, the values 
of Z(Y) are plotted as a function of Y In figure 10 for both Re and TaC. 

At Y = L, Z(Y) Is 10" 72 for Re and 10" 163 for TaC; W(Y) Is 10" 78 
g/cm 2 sec for Re and 10 -166 g/cm 2 sec for TaC. Thus, Insignificant 
vaporization occurs at the exit. Most of the vaporization takes place close 
to the entrance; here Z(Y) = 1, and W(Y) = W(nv). Therefore, the lifetime 
of a reslstojet tube should be governed by the lifetime of Its entrance; this 
Is the rate for normal vaporization. Our earlier calculations gave T"(Re) = 

42 hr and T"(TaC) = 5.6 sec. It Is also Interesting to compare the present 
findings with those for an evacuated tube. In an evacuated tube, the vapor- 
ization Is greatest at Its exit (fig. 9; eq. (52)) while for one with gas 
flow, the rate Is greatest at Its entrance (fig. 10; eq. (55)). 

For tubes with cross sections which are either nonuniform or noncircular, 
the values of A In equations (16) and (53) may be expressed In terms of Y. 
For most cases, the resultant equations are difficult to Integrate. However, 
equations (54) and (55) may always be used as an approximation. Regardless of 
the shape of the tube, the vaporization rate In a flowing gas-isothermal sys- 
tem does decrease as one proceeds from the entrance to the exit of the tube. 
This dependence of vaporization upon distance along the tube suggests a method 
for decreasing the actual vaporization of the tube. If a sacrificial sample 
of material Is placed at the entrance of the tube, the matrix gas could become 
saturated with the vaporized species, thus Inhibiting vaporization of the tube 
material, Itself. For TaC and other carbides, the sacrificial material could 
be carbon. 
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The average value of Z over the length of the tube can be obtained via 
equation (56) which was obtained by summing equation (55) over Y and then by 
dividing by L. 


Z(av) = ( sqr(M)/225 L) (1 - exp(-225 L/sqr(M))) 

= sqr(M)/225 L (56) 

The respective values of Z(av) for Re and TaC are 0.00607 and 0.00267. As 
expected, these are Identical to the values obtained earlier via equation (18). 

Stationary matrix gas . - When dependence on Y Is considered, one must 
start with the gaseous diffusion equation In a form similar to equation (28). 


F(Y) = -b M D(g) dP(Y)/dY (57) 

where F(Y) Is the flow of the vaporizing species In the Y direction. This 
value of F(Y) Is also related to P(Y) by equation (53). Combining these 
two equations results In a second-order differential equation. However, It Is 
readily Integrated and boundary conditions are easily assigned: the pressure 

of the volatile species at L Is taken to be zero. Thus, 


(P(v) - P(Y))/P(v) = 


cosh(Y sgr(Z sqr(2/birH)/R l D(q))) 
cosh(L sqr(Z sqr(2/birH)/R'0(g))) 


Via equations (3) arid (5) 


W(Y) = 


a Z P(v) sar(M) cosh(Y sor(Z SQr(2/birH)/R'D(q))) 
cosh(L sqr(Z sqr(2/birM)/R'D(g))) 


(58) 


(59) 


where Z again refers to the step coupled with the stationary gas step. By 
using equation (6), 


Z(Y) = Z cosh( Y sqr( Z sqr( 2/bwM)/R ' D(g) ) )/cosh( L sqr(Z sqr( 2/birM)/R ' D(g) ) ) 

= Z exp( ( Y - L) ( sqr( Z sqr(2/birH)/R , D(g))) 

= Z exp( 537 (Y - L) sqr(Z)) (60) 

The last form of equation (60) holds for TaC, If no chemical reaction Is In- 
volved. When the other step Is normal vaporization, Z = 1. For this case, 
calculated values of Z(Y) are plotted In figure 10. The vaporization behav- 
ior Is similar to that for an enclosure In that the rate Is greatest at the 
exit (Z(Y) = 1) and nil everywhere else. Thus, the lifetime of the tube 
equals the lifetime predicted for normal vaporization. At the entrance of the 
tube, Z(Re) = 10“ 154 ^ and Z(TaC) = 10 -233 ^. For nonuniform and noncircu- 
lar tubes, corrections could be made In equation (53) by expressing R' In as 
a function of Y. 


The average value of Z can be obtained by summing equation (60) over 
the length of the tube, L, and then dividing by L. The result Is 

Z(av) = 1/(2L sqr(Z sqr(2/birM)/R'D(g))) (61) 
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The values for Re and TaC are 1.4xl0 - ^ and 9.3x10“^, respectively. As 
should be the case, these are less than Z(max), 0.0025, calculated with the 
assumption that diffusion through the stationary matrix gas was the sole rate- 
controlling step. 

Dlffuslon/lnconoruent vaporization . - The coupling of vaporization with 
solid state diffusion controlled Incongruent vaporization requires special 
treatment. One begins with equation (33) but has to set different boundary 
conditions. We will use two different models. In Model 1, the vaporization 
rate of the carbon trlmer Is assumed to be proportional to the cube of the 
surface concentration of carbon, and the rate of vaporization of the carbon 
monomer to be proportional to the first power of Its concentration as In 
equation (47). 


W(ef ) = Z W(nv) = Z (C l3 W(C 3 ) + C' W(C)) (62) 

where Z Is the value for the step that Is being coupled with diffusion/ 
Incongruent vaporization. For normal vaporization Z = 1. This boundary con- 
dition makes Integration of equation (33) difficult so that the calculations 
were made via a variable finite difference method. The computer program Is 
given In table 1(b); the program need only be run until F(r) becomes less 
than 0.5. 

In Model 2, we assume that the rates of vaporization of both the carbon 
trlmer and the carbon monomer are proportional to the first power of the car- 
bon surface concentration. This boundary condition allows easy solution of 
equation (33) (ref. 7, p. 56). For planar diffusion and with our symbolism. 


w 

F(r) = 2E 2 £ exp(-B(n) 2 D(s) T'/S)/B(n) 2 (B(n) 2 + E 2 + E) 
n=0 


(63) 


and 


Z(ef ) - C' - 2E £ ( -B( n) 2 D(s) T'/S)/(B(n) 2 + E 2 ♦ E) 
n=0 


(64) 


where 


B(n) tan B(n) = E 


(65) 


and 


E = Z E(nv) = Z S W(nv)/D D(s) = a Z S P(v) sqr(M)/D D(s) (66) 

The Z again refers to the value for the step being coupled with diffusion/ 
Incongruent vaporization. When the step Is normal vaporization, Z = 1. The 
computer program for Model 2 Is given In table 1(c). 

The values of F(r) for Models 1 and 2 are plotted In figure 6 as a 
function of the square root of time. The curves Indicate roughly parabolic 
dependence. At all times, the values of F(r) are greater than those 
predicted by Model 0. This Is In accord with our earlier comment that the 
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effective rate of two combined series steps Is always less than either one 
alone. From figure 6, the value of T" for Model 2 Is 179 sec, not much more 
than that for Model 0 (T“ = 162 sec). However, T" for Model 1 Is 309 sec, 
definitely greater. Thus, one sees that the value of T" Increases with the 
order of the vaporization step; remember that the order for Model 0 Is zero 
Plots of Z(ef) In figure 7 confirm the fact that the vaporization rate for 
Model 0 Is the greatest, the rate for Model 2 Is slightly less, and the rate 
for Model 1 Is the least. For Models 1 and 2, carbon concentration profiles 
(not Included In this report) were found to be similar to those for Model 0 
(fig. 8). Moreover, the half lives for Models 1 and 2 were both found to be 
proportional to the square of the thickness, S, of the reslstojet tube. One 
can again calculate via equation (36) a value of Z' (a pseudo values for 
Z). The value of Z' for Model 1 Is 0.018; that for Model 2 Is 0.0031. 

Chemical vaporization . - Equation (38) can be used to give the value of 
Z(ef) for the coupling of normal vaporization with chemical vaporization. 
These are parallel steps. The resultant value of Z(total) for TaC Is, 
therefore, 578. No chemical reaction Is presumed for Re; therefore, Z(total) 
would be unity. 


Coupling Chemical Vaporization With Other Steps 

Coupling I ndependent of Y and T 1 . - Because Z(cs) » 1, coupling of 
chemical vaporization Is similar to the coupling of normal vaporization. In 
other words, the rate of chemical vaporization Is slowed down by the other 
steps. This In turn allows the chemical reaction Involved to proceed closer 
to equilibrium, and the equations used will therefore have greater validity. 
To avoid the difficulty of defining Z values In equation (48), we begin as 
follows: 


Z(ef) = W(ef)/W(nv) = W(ss)/W(nv) (67) 

Remember that "ss" refers to the slowest step. In the calculation of W(ss), 
one uses the P, M, and M' for chemically formed species rather than for the 
normally formed species. Thus, for a stationary matrix gas 

Z(ef) = Z(ss) P(cs) M'(cs)/P(nv) sqr(M(nv) M(cs)) 

= Z( ss) Z(cs) ( 68 ) 

Remember Z(cs) = 577. For a flowing matrix gas and the boundary layer ef- 
fect, equation (67) leads to 

Z(ef ) = Z(ss) P(cs) M(cs)/P(nv) M(nv) (69) 

For TaC, this Is approximately the same as equation (68). However for the 
coupling of coating, Z(ef) = Z(ss) because diffusion through a coating Is 
Independent of the values of P and M of the substrate material. Thus, 

Z(cs + flowing gas) = 1.54; Z(cs + boundary layer) = 16.1; Z(cs + stationary 
gas) = 1.44 and Z(cs + coating) = 1.43xl0‘ 6 . The values of Z for the 
flowing gas and the stationary gas are expected to be functions of Y. Thus, 
the values just given are only average values. 
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Flowing gas . - To determine Z(Y) as a function of Y, the first form of 
equation (55) can be used. The Z for C2H2 In this equation Is taken as 
577 (Z(cs) ) and M as 24.02 (M' (C 2 ) ) - The values of Z(Y) from this 
equation are plotted In figure 10 as a function of Y. One sees that with 
chemical vaporization as well as with normal vaporization, the greatest vapor- 
ization rate occurs near the entrance of the tube. At Its exit, the vaporiza- 
tion rate Is nil while at the entrance, the rate Is equal to the rate for 
chemical vaporization alone (Z(Y) = 577). The value of. Z(Y) for CH4 Is 
small, 2.8, with respect to the Z(Y) for C2H2 at Y = 0, but It falls 
off less rapidly with Y (eq. (55)). It Is Interesting to note that except 
for small Y values, the rate of vaporization of C3 exceeds that of 
C2H2 because the matrix gas becomes more quickly saturated with C2H2 
than with C3. Again from the practical standpoint, the lifetime Is still 
controlled by the lifetime at the entrance, and the dependence of Z(ef) on 
Y Is of secondary Importance. 

As In the case of normal vaporization, corrections In the equation can be 
made for noncircular and nonuniform diameter tubes. Also, chemical vaporiza- 
tion could be Inhibited by placing a sacrificial vaporizing material such as 
carbon at the entrance to the tube. This would saturate the matrix gas with 
C2H2 and prevent the chemical vaporization of the tube Itself. 

Stationary matrix gas . - Via equation (60), one finds the dependence of 
Z(Y) on Y. The Z In this equation Is Z(cs), 577; the M Is M'(C2), 

24.02 g/mole. The calculated values of Z(Y) are plotted In figure 10. At 
the exit of the tube Z(Y) Is 577. It Is nil everywhere else, being lO -2 ^ 
at the entrance. The life of the tube Is controlled by the life of the exit. 

Dlffuslon/lncongruent vaporization . - To couple chemical vaporization 
with this step, we return to Models 1 and 2. The corresponding models that 
Include this chemical vaporization contribution are referred to, respectively, 
as Models XI and X2. The boundary condition for Model XI Is obtained by modi- 
fying equation (62) to read 

w = c' 2 w(c 2 h 2 ) + c w(ch 4 ) + c' 3 w(c 3 ) + r w(C) (70) 

In equation (70), the first two terms give the chemical contribution while the 
last two which are relatively Insignificant give the normal vaporization con- 
tribution. For Model X2, equation (64) may be used by setting value of Z In 
equation (66) equal to 577 and M equal to 24.02. 

The values of F(r) and Z(ef) for Models XI and X2 can be calculated 
via the programs In table I. The values are plotted respectively In figures 6 
and 7. The curves correspond very closely to the respective curves for Model 
0. This means that the vaporizable species Is vaporizing as soon as It Is 
formed on the surface; It Is due to the high vapor pressure of the chemically 
formed species. Rapid vaporization Is the chief postulate of Model 0. The 
fact that the chemical vaporization rate does not exceed the rate predicted by 
Model 0 emphasizes that Model 0 predicts the fastest rate possible for the 
process If controlled by solid state diffusion. 
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Other Interactions with the Dlffuslon/lncongruent Vaporization Step 

Coupling Independent of Y . - Again we begin with Models 1 and 2. To 
Incorporate other rate-controlling steps Into these models, we again modify 
the boundary conditions. The boundary condition for Model 1 Involves the def- 
inition of W. Equation (62) gives the desired condition If Z is the Z 
value for the particular step being coupled with dlffuslon/lncongruent vapor- 
ization. For models derived from Model 2, the value of E (eq. (66)) deter- 
mines the boundary conditions. The Z In this equation Is the Z for the 
step being coupled with dlffuslon/lncongruent vaporization. To distinguish 
these models from the original one, we refer to the ones Including the bound- 
ary layer contribution as Models B1 and B2. Those Including the surface coat- 
ing, we call Cl and C2; those Including the evacuated enclosure, El and E2; 
and those for a flow gas, FI and F2. Because the value of Z for the sta- 
tionary gas Is the same as the Z for the evacuated enclosure. Models El and 
E2 will be used to represent the stationary gas effect. 

The basic computer programs In table I can again be used to calculate the 
value of F(r) and Z. These calculated values are plotted as a function of 
time, T' t In figures 11 and 12, respectively. In all the cases at a specified 
time, the carbon concentration F(r) for the coupling of rates Is greater 
than that predicted for normal vaporization (fig. 6). Eventually, the carbon 
concentration for all models will go to zero. From figure 11, we can obtain 
the values of T" for the different models. The values of T" are as fol- 
lows: Model B1 (boundary layer). 1200 sec (20 min); Model B2, 470 sec (7.8 
min); Model Cl (coating), 1.2xl0 7 sec (139 days); Model C2 (coating), 

5.1x10° sec (59 days); Model El (evac. enclosure or sta. gas), 7.6xl0 3 sec 
(2.1 hr); E2 (evac. enclosure or sta. gas), 3.3xl0 3 sec (0.92 hr); Model FI 
(flowing gas), 7.1xl0 3 sec (2.0 hr); and Model F2 (flowing gas), 3.1xl0 3 
sec (0.86 hr). In all these cases, the value of T" for the models (Bl, Cl, 
etc.) Involving third-order dependence on concentration Is over 2 times larger 
than that for the corresponding models (B2, C2, etc.) Involving first-order 
kinetic dependence. The values of Z via equation (36) are: 0.0047 for 81, 

0.0119 for B2 , 4.7xl0~ 7 for Cl, l.lxlO" 6 for C2, 7.4xl0” 4 for El, 

1.7xl0- 3 for E2, 7 .9x1 0 -4 for FI, and 1.8xl0” 3 for F2. 

In three of these coupling conditions (evacuated enclosure, flowing gas, 
and stationary gas) the values of Z(ef) are expected to be a function of 
Y. Thus values of Z(ef) obtained via Models El, E2, FI, and F2 are only av- 
erage values. For the flowing gas, the value of Z at the entrance Is great- 
er than the average Z from either Model FI or F2. And for the stationary 
gas and evacuated enclosure, the Z at the exit Is greater than the average 
Z from either Model El or E2. We will discuss the dependence of Z on Y In 
the next section. 

Evacuated enclosure . - The coupling with the evacuated enclosure Is a 
relatively simple problem to treat as a function of Y. We shall call the 
Models E'l and E'2. The boundary condition for Model E'l Is obtained by set- 
ting Z In equation (52) equal to one, and substituting the value of Z(Y) 
from equation (52) for the value of Z In the boundary condition (eq. (62)). 
This results In a boundary condition equation which Is dependent on Y. 

W = W(nv) (1 - cos(arctan(R'/(L - Y))))/2 (71) 
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The boundary condition for Model E'2 Is obtained by setting the Z In equa- 
tion (51) equal to unity and substituting the value of Z(ef) for the value 
of Z In the boundary condition In equation (66). 

E = E(nv) (1 - cos(arctan(R'/(L - Y))))/2 (72) 

The calculations via the programs In table I yield F(r,Y), the average carbon 
concentration at a position Y along the tube (enclosure). Calculated values 
of F(r,Y) and Z for Model E'2 are plotted as a function of T‘ In figures 
13 and 14, respectively. Me have also plotted the result from Models El and 
E2 In these figures for comparison. The values obtained for Model E2 corres- 
pond roughly to a Y value of 9. To prevent a cluttered graph, the values 
for Model E'l are not plotted. The values for Model E'l vary from those for 
Model E'2 In roughly the same manner as those for Model E2 vary from Model 
El. From a more detailed graph than that given In figure 13, we obtained the 
values of T" as a function of Y and have plotted them In figure 15. This 
figure gives a vivid picture of how the relative lifetime varies along the 
tube. Except for portions of the tube near the exit, values of T" are rela- 
tively large, but the life of a tube must be judged by shortest lifetime of 
any of Its parts. 

Flowing matrix gas . - The Y dependence of the coupling of the diffusion/ 
Incongruent vaporization step with the existence of a flowing matrix gas is 
relatively complicated. A two dimensional computer program would result which 
would be difficult to run In a reasonable time on a small computer. Moreover, 
the two dimensional results would be difficult to Illustrate In graphic form. 

In addition, the calculation In most cases would be only of academic Interest 
because all we should be Interested In Is the lifetime of the entrance which 
has already been shown to be the position of shortest lifetime. The curves 
for Models 0, 1, and 2 In figure 7 can be used to Illustrate the vaporization 
at the entrance as a function of time. 

Instead of setting up a rigorous model to determine the dependence of Z 
on Y, we set up a crude model. We will not bother the reader with the calcu- 
lations or plots. This model shows that at the beginning, the greatest vapor- 
ization rate exists at the entrance. But as time proceeds, this region Is 
depleted of carbon. Thus, the position with the greatest vaporization rate 
proceeds with time along the tube. The position of the maximum rate Is rough- 
ly proportional to the time. Except for very short times the vaporization 
profile relative to the position of maximum vaporization does not change much 
with time. In other words, the vaporization profile proceeds along the tube 
In a slow single pulse fashion. 

Stationary matrix gas . - A similar result would be expected for the Y 
dependence with a stationary matrix gas, except the maximum rate would begin 
at the exit of the tube and proceed toward the entrance with time In a slow 
single pulse manner. The dependence on time of Z at the exit for the sta- 
tionary gas can be depicted by the curves for Models 0, 1, and 2 In figure 7. 


Coupling of Surface Coating 

Coupling Independent of Y . - Four steps that can couple with the coating 
(ct) have not yet been discussed: gaseous boundary layer, evacuated enclo- 
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sure, flow matrix gas, and stationary matrix gas. In all these cases the 
value of Z Is much greater than the Z for the surface coating so that the 
value of Z(ef) = Z(ct). Therefore equation (32) Is directly applicable; for 
Re, Z(ef ) = 1. 6xl0- 3 , and for TaC, Z(ef) = 1.43xl0" 6 . Note that Z(ct) 

Is based on many hypotheses and an assumed value for the thickness of the 
coating. Also, remember that no coating Is contemplated for the reslstojet. 

Coupling dependent on Y . - The equations that can be used have already 
been presented. For coupling with the enclosure one uses equation (52); for 
coupling with a flowing gas, equation (55); and for a stationary gas, equation 
(60). In all these equations, the Z Is the one for the surface coating. 

The calculated values of Z(Y,Re) and Z(Y,TaC) for the enclosure are plot- 
ted In figure 9, those for the flowing gas are given In figure 16, and those 

for the stationary gas are given In figure 17. 

Because the Z for the surface coating effect Is quite small, the varia- 
tion of Z(Y) with Y Is small In some cases, particularly for the station- 
ary matrix gas. Nevertheless, for the stationary gas as well as for the evac- 
uated enclosure, the value of Z(Y) Is always somewhat greater at the exit 
than at the entrance, and the Z(Y) for the flowing gas Is always somewhat 
greater at the entrance than at the exit. However, this variation for Z(TaC) 
cannot be seen In figures 16 and 17 due to the scale of the graph. 

In all these cases, we are Interested primarily In the maximum value of 

Z(Y) which we have just discussed; the maximum rate determines the life of 

the entire tube. It Is at the entrance for the flowing gas and Is at the exit 
for the stationary gas where Z(Y) = Z(ct). The maximum value at the exit for 
the evacuated enclosure; Z(Y) = Z(ct)/2. From Part I, Z(ct,Re) = 1.43xl0 -3 
and Z(ct,TaC) = 1.6xl0~&. If a coating exists, the life of the tube Is 
always controlled by the diffusion through It and not by the steps Involving 
an evacuate enclosure, a stationary gas, or a flowing gas. 


Coupling Involving the Boundary Layer 

Couplings Independent of Y . - In the case of coupling the boundary layer 
effect with a flowing gas and a stationary gas, we find that the Z for the 
boundary layer Is larger than the other two. If a much larger thickness were 
assumed for the boundary layer, this would not be true. However, with our 
assumptions, Z(av) Is equal to Z for the process being coupled with the 
boundary layer effect. For the coupling with the flowing gas, Z(av,Re) = 
0.00607, and Z(av,TaC) = 0.00267. For coupling with the stationary gas, 
Z(max) = 0.0025 for both Re and TaC. 

Coupling dependent on Y . - Again, we can use equations (55) and (60) 
respectively for the coupling of the boundary layer effect with the flowing 
matrix gas and with the stationary matrix gas. The Z used Is the Z for 
the boundary layer (Bl) diffusion. The calculated values for Re and TaC are 
plotted respectively In figures 16 and 17. As In the other examples of 
coupling, Z(Y) for the flowing gas Is greatest at the entrance; It Is equal 
to Z(bl). The Z(Y) for the stationary gas Is greatest at the exit where It 
Is also equal to Z(bl). Again, the life of the tube In a flowing gas Is de- 
termined by the life of Its entrance and the life In a stationary gas Is 
determined by the life of the tube's exit. In both cases, the value of 
Z(max) Is set by the value of Z(bl). 
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PART III: PRACTICAL CONSIDERATIONS IN ISOTHERMAL SYSTEMS 

In Part IV, we will see that the existence of a thermal gradient In a 
practical system can tremendously Influence the vaporization rate of the ma- 
terials Involved. Before we get Into the variations that a temperature gradi- 
ent gives rise to, we desire to discuss the practical aspects of a purely Iso- 
thermal system because some practical systems may have to be considered as 
Isothermal . 

Part II has shown us that for Isothermal systems, equations can be de- 
rived for the case of coupled kinetic steps. We limited our discussion of 
such equation to two-step processes to prevent verboseness In the presentation. 
As we will see In the following discussion, coupling two step is generally 
sufficient for a particular problem. Coupling of all the steps presented Is 
never needed, and some combinations of steps cannot be coupled together. For 
example, a flowing gas effect cannot be coupled to an evacuated enclosure ef- 
fect. Now that we have examined the equations for the single and double step 
processes, we are In a position In Part III to see which steps are really Im- 
portant for an Isothermal problem. 


A Look at the Kinetic Steps 

Enclosures . - For this discussion, one can group together all the steps 
that are concerned with an enclosure; these Include the evacuated enclosure, 
and an enclosure with either a flowing or stationary gas. As we have seen, 
each of these steps add a dependence upon Y. However, the life of the con- 
tainer (tube) Is determined by the shortest lifetime of any part of tube. 
Moreover, this shortest lifetime Is not dependent on the existence or dimen- 
sions of the tube or other enclosure. Thus, the existence of an enclosure 
whether evacuated or containing a matrix gas has essentially no effect on the 
effective lifetime of the tube. One has to consider only the position of the 
maximum rate of vaporization which Is either at the entrance or exit depending 
upon the specific conditions. Remember, we are discussing an Isothermal sys- 
tem. The existence of a temperature gradient greatly affects the calculations 
In systems Involving enclosures. 

Vaporization per se . - We may make another grouping of kinetic steps. 
Normal vaporization and chemical vaporization fall Into a common category. 

They are parallel steps and thus do not Interact with each other. Also, chem- 
ical vaporization Interacts with every other step In the same manner as the 
normal vaporization step. The effective rate of their combination Is merely 
the sum of their Individual rates. Moreover, one of the two steps Is general- 
ly more Important than the other so that one can generally Ignore one of these 
two parallel steps. For Re, chemical vaporization can be Ignored. For TaC, 
normal vaporization can be Ignored. But In any vaporization problem, one of 
these two processes must be considered. Even when a coating makes the vapor- 
ization step of the substrate unimportant, the vaporization of the coating 
Itself must be considered as discussed In the next paragraph. 

Coating . - For the reslstojet, no coating Is contemplated. In order for 
a coating to Increase the life of the system, the coating would have to have a 
vapor pressure lower than that of the substrate. However, If It had such a 
low pressure, then the tube should probably be made from the coating material. 
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Even If a coating were used In a particular situation, the Important vaporiza- 
tion problem would probably Involve the vaporization of the coating and not 
the loss of substrate via diffusion through the coating. Nevertheless, a 
coating can certainly be used to shield from vaporization a very volatile ma- 
terial that has to be present In the system. Remember "that our earlier calcu- 
lations show that a coating can be quite effective In Inhibiting the 
vaporization rate of a substrate. 

Olffuslon/lncongruent vaporization . - If the material under consideration 
vaporizes Incongruently, then this would always be an Important step to con- 
sider In the overall process. If the vaporization Is congruent (e.g., Re), 
consideration of this step does not arise. 

Boundary layer . - In any system with a flowing matrix gas, a boundary 
layer should exist and have the effect of Inhibiting the rate of vaporization 
as seen earlier. Thus, this effect should Indeed be considered In the calcu- 
lation of vaporization rates for the reslstojet concept. Only for concepts 
Involving vaporization Into vacuum or very low pressure, should consideration 
of the boundary layer be Ignored. 


Application to the Reslstojet 

Use of Re . - By reflection on the above discussion, one can see that the 
steps we should consider are normal vaporization and diffusion through the 
boundary layer (bl). The value of Z for this combination has already been 
determined; Z(bl) = 0.0607 for a 0.005 cm thick boundary layer. Thus, W = 
4.2x1 0“° g cm 2 /sec at the entrance of the tube. More accurate calcula- 
tions can be made If a more accurate value of the boundary layer thickness Is 
known. At 665 hr (28 days, value of T") for a 0.01 cm thick tube at tempera- 
ture, we predict that half the thickness of the tube at the entrance will be 
evaporated. The life of Re can be extended still more by the use of the sac- 
rlflcal material discussed earlier In this report. 

In a practical problem, one must also look at more than the time during 
which the tube Is full of gas and Is at the temperature of operation. Of 
course, when the tube Is at room temperature the vaporization rate is nil. 
However, during the heating and cooling time, the tube may not always contain 
a matrix gas. The rate of vaporization for such periods Is determined by the 
normal vaporization rate of the material which at temperature Is faster by a 
factor of 1/Z(bl) than the rate for when the boundary layer of a matrix gas Is 
present. As long as the time at temperature with no gas Is small with respect 
to the time with the matrix gas, one need not be concerned with this situation 
for the reslstojet. Also, to our advantage Is the fact that without the ma- 
trix gas, the most rapid vaporization rate occurs at a different part of the 
tube, the exit. 

As far as vaporization Is concerned. Re appears to be a very good candi- 
date for the exit tube of the reslstojet and other extremely high temperature 
devices, especially In situations where chemical reactivity can be Ignored. 

Use of TaC. - For an evaluation of the vaporization behavior of TaC, one 
should consider the vaporization mode Itself to be chiefly chemical vaporiza- 
tion In nature. Incongruent vaporization Is also an Important step as well as 
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diffusion through the gaseous boundary layer. In the treatment of this prob- 
lem, we will restrict ourselves to a modification of Model 1, which we believe 
Is closer to reality than a variation of Model 2. We will refer to the model 
as Model X'l. The boundary condition Is a variation of equation (70). 

W(ef ) = Z(C' 2 W(C 2 H 2 ) + C' W(CH 4 ) + C' 3 W(C 3 ) + C' W(C)) (73) 

The value of Z would be that for the boundary layer effect. The program In 
table 1(b) can again be used. 

The calculated values are shown In figure 18. From this figure, the 
half-life, T", of TaC Is 192 sec under these conditions. This Is the time 
required for the average composition of the TaC at the entrance to decrease to 
Ta 2 C. This value of 192 sec Is not long enough for TaC to be considered as 
a candidate for the exit tube of the reslstojet. Moreover, the surface con- 
centration at the entrance reaches 50 percent depletion In about 2 ys (fig. 

18) and over 97 percent depletion at 192 sec, the value of T". At these de- 
grees of depletion, the chemical and physical properties of the material are 
no longer that of TaC. Thus, the value of T" does not represent the life- 
time of the physical properties. Such a lifetime may be much less. 

There are still other factors to consider. Many of the comments made for 
the Re also holds for the TaC. Thus, the rate of vaporization Is different 
during the heating and cooling periods when no matrix gas Is present. How- 
ever, for TaC at temperature with no matrix gas, the net rate of vaporization 
per time Increment will be less because rapid chemical reaction will not occur. 
Because we may also apply the principle of sacrificial material In the case of 
TaC, the rate of vaporization of the TaC tube. Itself, could be decreased con- 
siderably as we have mentioned several times before. Thus, the possible use 
of TaC In the reslstojet depends upon whether and how much sacrificial materi- 
al such as carbon could be Incorporated Into the design of the reslstojet. 

One must also remember that In our calculations we have assumed that the 
chemical reactions went to completion. We have assumed the worst situation. 
Under the conditions of the reslstojet, chemical vaporization might not occur 
at all. This can only be determined by experiment. It has been reported 
(ref. 10) that under a static hydrogen at a half atmosphere, no reaction be- 
tween TaC and H 2 could be detected. Even If no reaction occurs the value 
of T" (Model B1 ) would only be 20 min, far less than the value for Re. 


PART IV: EFFECT OF A TEMPERATURE GRADIENT 

Importance of a Temperature Gradient 

In Part IV we shall see that the existence of a temperature gradient In a 
system may be one of the most Important factors that should be considered for 
some vaporizing systems, particularly that of the reslstojet. Because the 
existence a temperature gradient does not constitute a kinetic step, discus- 
sion up to this time has been concerned with only a single temperature. For 
single kinetic steps, the existence of a temperature gradient would merely 
mean that the rate of vaporization would be different at different locations. 
This Is particularly evident If we consider normal vaporization. The vapor- 
ization rate would vary In an exponential manner with the absolute temperature, 
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but vaporization at one point In space would have no effect upon the vaporiza- 
tion at another point In space. 

However, In some cases Involving coupling of kinetic steps, the vaporiza- 
tion rate at one location In space can have an effect upon the rate at another 
location. This arises from the Interaction of two relationships: (1) the 

vaporization rate Is proportional to the difference between the vapor pressure, 
P(v), and the partial pressure, P(Y), of the vaporized species In the matrix 
gas; and (2) the value of P(Y) can be affected by the value of P(v) else- 
where In the system where the temperature can be different. (If P(Y) were 
greater than P(v) In some cooler location, condensation would even occur). 
However, the Interaction can only occur In some step Involving an enclosed 
system. In an evacuated enclosure, the Interaction does not occur because 
P(Y) = 0. Moreover for the case of a stationary matrix gas In an enclosure, a 
temperature gradient should not affect the vaporization at the exit where 
P(Y) = 0. And It Is here where the vaporization rate Is the greatest, and 
thus this Is the place that determines the lifetime of the entire tube. 

Therefore, when evaluating the effect of a thermal gradient, we need con- 
sider only the coupling of the flowing matrix gas In an enclosure with other 
relevant kinetic steps. This, of course, Is the situation existing In the 
reslstojet and In many other practical systems. 


Introduction of the Temperature Gradient Effect 

Many of the parameters used In this report are dependent on the tempera- 
ture of the material. The existence of the coefficient of expansion makes one 
realize that even the radius of the tube Is a function of the temperature. 

The use of all of the temperature dependent terms would become cumbersome and 
quite Involved. Luckily, most of the dependencies on temperature can be Ig- 
nored because they are Insignificant compared with the dependence of the vapor 
pressure on temperature; this Is exponential In nature. The vapor pressure of 
a substance Is theoretically related to the absolute temperature as follows: 

P(v) = k-, exp (-k 2 /T) (74) 

where the k's are constants for the element or compound under consideration. 
To permit easier mathematical treatment we have approximated equation (74) as 
follows: 


P(v) = k-, exp(k 3 T) ( 75 ) 

This Is permitted because we are Interested In only relatively small differ- 
ences In temperature. Because the temperature dependence of the vapor pres- 
sure of Re and the Co component of TaC are similar (refer to the slopes of 
the curves In fig. 2), k 3 In equation (75) can be approximated to be the 
same value for both, 0.012 reciprocal degrees. If a uniform temperature gra- 
dient, G', Is assumed to exist along the reslstojet tube, the vapor pressure 
can easily be expressed as function of Y. 

P(v,Y) = P(v,2773) exp(0.0l2 G' (Y - L)) (76) 
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where 6' Is In units of degrees/cm, and the equation assumes the temperature 
at the exit, L, Is 2273 K. 

As discussed In the first section of Part IV, the effect of a temperature 
gradient In most situations Is merely to produce a different rate of vaporiza- 
tion at the various parts of the tube. The calculations of such effects are 
straightforward, and thus no further detail should be required. However, when 
the rate a vaporization Is dependent on P(Y), the applicable equations are 
modified considerably. In the following sections, we Illustrate the modifica- 
tion of the kinetic equations for the practical situation of a flowing matrix 
gas. 


Coupling With a Flowing Gas 

The first quantity that Is desired Is, of course, the partial pressure of 
the vaporized gas as a function of Y because this Is the quantity that In- 
hibits the vaporization rate. This can be obtained by combining equations 
(3), (16), (53), and (76). If P(Y) =0 at Y = 0, then, 

pm exp(k 3 G'Y) - exp(-225 Z Y/sqr(M)) 

P(v,2773) = ex P(- k 3 6 ' L ) (1 + k 3 G' sqr(M)/225 Z) (77) 

where 225 Is R' sqr(2ir/b)/F' . If P(Y) equals the vapor pressure of the 
material at position Y = 0, that Is 

P(Y) = P( v,2773) exp(k 3 G'L) (78) 

then 

pm exp(k G'Y) + k G'sqr(M) exp(-225 Z Y/sqr(M))/225 Z 

— T ' - pvn( _k fi'll 2 

P( v,2773) exfn k 3 (1 + k 3 G'sqr(H)/225 Z) 


(79) 


The Z In these equations Is the value for the rate-controlling step. For 
normal vaporization, Z = 1; for a boundary layer controlled process In a re- 
slstojet with our parameters, Z = 0.028 for TaC. For both equations (77) and 
(79), the relative rate Z(Y) Is given as follows: 


Z(Y) = Z (P(v,Y) - P( Y) )/P( v,2773) 

If equation (77) Is assumed to be valid, then 

zm - Z exof 0 012 g 1 LI H ex P (°- 012 G'Y) + exp(-225 Z Y/sqr(M)) 
17 - ^ exp(-u.ui<: t. L) (1 + o.012 G'sqr(M)/225 Z) 


where H = 0.012 G'sqr(H)/225 Z. 

However, If equation (79) Is a better approximation, then 


Z(Y) = Z exp(-0.012 G'L) 


exp(0.012 G'Y) - exp( -225 Z Y/sar(H11 
(1 + 225 Z/0.012 G'sqr(M)) 


(80) 


(81) 


(82) 
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The values of Z(Y) for Re and TaC from equation (82) are plotted as a func- 
tion of Y In figure 19 as dashed lines. Normal vaporization Is assumed to 
be the rate-controlling step, and the value of G' Is taken as lOVcm. The 
solid line segments Indicate how Z(Y) In equation (81) varies from Z(Y) In 
equation (82). The difference In the two equations Is due to the postulate 
concerning the pressure of the vaporizing species at Y = 0. The curves for 
equation (81) show that the rate of vaporization Is greatest at Y = 0; this 
Is a reflection of the postulate that the pressure of the volatile species at 
Y = 0 Is zero. The value of Z(Y) at Y = 0 for equation (82) Is zero but 
could not be depicted on figure 20 because the abscissa Involves a logarithmic 
scale; the zero value Is a reflection of the postulate that P(Y) at Y = 0 
is taken to be the vapor pressure at that position. 

In a system like the reslstojet, the tube leading Into the exit tube would 
probably be made out of the same material as the exit tube Itself. At least 
It would be a desirable situation because In such a situation the matrix gas 
could contain the vaporizing species near the saturation value for the entrance 
temperature of the exit tube. Under this condition, there would be very little 
vaporization occurring at Y = 0. But this Is difficult to express mathemati- 
cally. Nevertheless, the solution for the desired assumption should be the 
straight dashed lines In figure 19 with values of Z(Y) at Y = 0 being on 
the respective lines. Note that In the nonlsothermal model, the greatest 
vaporization rate occurs at the exit, the position of highest temperature, 
rather than at the entrance as Is the case for the Isothermal model. 

From a practical viewpoint, figure 19 Is one of the most Informative 
graphs In this report. It tells us that the vaporization rate at any spot Is 
not as severe as predicted directly from an Isothermal rate-controlling equa- 
tion. The fact that the gas space already contains some vapor of the vapor- 
izing species Inhibits the rate. The vapor already present arises from the 
vaporization of the material In cooler portions of the system, that Is, It Is 
due to the presence of a temperature gradient. For Re, the average rate Is 
about 1/200 of that predicted by normal vaporization; for TaC, the rate Is 
about 1/500 of that predicted by normal vaporization. Remember these values 
are for a gradient of 10°/cm. 

The next step Is to show the effect of G 1 on the profile. Figure 20 Is 
based on equation (82) and normal vaporization of TaC. For a very small gra- 
dient, G 1 = 0.1°/cm, the average rate Is very low because the saturation value 
at Y = 0 is close to the saturation at Y = L. Also, the curve Is almost 
horizontal, that Is, the rate Is relatively uniform from entrance to exit. 

For large values of gradient, G' = 100 o /cm, the profile Is not flat; this 
arises from the exponential dependence of Z(Y) on the temperature and there- 
fore on Y. However even at Y = L, the value of Z(Y) does not approach 
unity. Moreover, the existence of a large value for G' Is not expected In 
the reslstojet. 

By using the proper values for the parameters In equation (82), the val- 
ues of Z(Y) for various rate-controlling steps can be ascertained. Figure 
21 shows the effect of flow rate as well as the effect of boundary layer dif- 
fusion being the rate-controlling step. From this figure, one can Infer that 
Z(Y) Is Inversely proportional to the flow rate of the gas In the range under 
consideration. Also, one can Infer that the value of Z(Y) for boundary lay- 
er diffusion Is directly proportional to Its value of Z; Z(bl) = 0.028 for 
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TaC. Thus, the various kinetic steps which affect the vaporization rate for 
an Isothermal process can have the corresponding effect on vaporization when a 
temperature gradient Is present. 

From the boundary layer curve In figure 21, It Is not unreasonable to 
expect that materials for the reslstojet could last more than 1000 times long- 
er than that predicted by simple Isothermal vaporization. For TaC, the exis- 
tence of a chemical reaction can negate much of the effect caused by the 
thermal gradient: Z(Y,cs) = Z(Y,nv) x Z(cs). Thus the Z(Y)'s In figures 21 

and 22 should all be multiplied by 577 to obtain Z(Y,cs). 

The effect of the temperature gradient on Incongruent vaporization may 
also be considered. We will call It Model X"1 . One starts with Model 1 (see 
eq. (33)) using equation (73) as the boundary condition. Equation (73) In- 
cludes the contribution from chemical vaporization. The Z In equation (73) 
would be the maximum Z(Y) value from equation (82), that Is the value at Y 
= L. If the boundary layer Is considered to control the rate, then Z In 
equation (73) would be 8.05xl0 - 5, obtainable from the boundary layer curve 
In figure 21. Via a modification of the computer program In table 1(b), the 
values of F(r), Z, and C' were calculated as a function of time and plotted 
as a function of time In figure 22. The value of T" Is found to be 660 sec; 
this Is not sufficiently long to expect TaC to be a good material for the re- 
slstojet exit tube. ■ 

However, Re was found to be a good prospect for the reslstojet tube under 
Isothermal conditions. It would have a lifetime of 665 hr. It Is now ex- 
pected to last even longer under thermal gradient conditions, perhaps over 
100 000 hr. Of course, this presumes that no unforseen process Is operative, 
such as chemical vaporization. 


Another Approach 

If the value of G' Is small, we may be satisfied with an average value 
for Z. In Part I of this report, equation (18) (flowing gas equation) was 
derived by using the assumption that the gas leaving the tube was saturated 
with the vaporizing species, but that the gas entering contained no vaporizing 
species. To correct this equation for the fact that the gas entering could be 
saturated with the species at the temperature of the entrance, we sho uld mul- 
tiply by the factor (P(v,2773) - P(v,2773-L G' ) )/P( v,2773) . This factor Is 
the ratio of the pickup In the tube divided by the total pickup In the system. 
Through the use of equation (76), this factor can be shown to be equal to 1 - 

exp(0.012 G'L). Thus, for our reslstojet parameters. 

r 

Z(av) = (1 - exp(0.012 G'L))(b F'sqr(M)/a A(s)) 

= 4.45X10- 5 (1 - exp(0.1 2 G')) (83) 

For small values of G' 

Z(av) = 5.34x10-5 G'sqr(M) (84) 

For TaC with a G' value of 0.1°/cm, Z(av) = 3.2x10-5; this Is essentially 
the value obtained previously via equation (82) by using the same G‘. Refer 
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to the lowest curve In figure 20. To find Z(av) for a boundary layer dif- 
fusion controlled process, It Is logical to multiply this Z(av) by the value 
of Z for the rate-controlling step (0.028 for the boundary layer with TaC). 

Our first approach to the temperature gradient effect Is more stringent 
than this latter one. However, this latter approach Is simpler and Is more 
easily applied than the former. Both methods point out that the existence of 
a thermal gradient along a tube with a flowing gas can spread out the vapor- 
ization effect over the entire area of the tube. Thereby, It can prolong the 
life of the tube as compared with the life predicted from simple Isothermal 
calculation. Thus, many materials will perform much better (by perhaps a fac- 
tor of 1000) than predicted by simple vaporization theory. 


4 SYNOPSIS 

In Part I, we discussed the Individual effects of various kinetic steps. 
In some cases chemical reaction can severely enhance the rate. The existence 
of a gas In an enclosure slightly Inhibits the rate. In some cases, the use 
of a diffusion barrier may be advantageous. 

In Part II, the coupling of kinetic steps Is considered. Parallel steps 
Increase the rate; series steps decrease the rate. For some situations, coup- 
ling steps also Introduce a dependence on position. Dependence on position 
suggests the use of a sacrificial vaporizing material. 

In Part III, vaporization Is viewed In practical Isothermal systems. 

When a matrix gas Is Involved, the existence of a boundary layer Is a very 
Important factor to consider. Other Important factors to consider (such as 
chemical reaction) depend upon the materials under consideration. Rhenium 
appears to be a likely candidate for use at 2500° C In the reslstojet but tan- 
talum carbide does not. 

In Part IV, the existence of a temperature gradient Is considered. This 
Is found to be a very Important factor. The material vaporized In low temper- 
ature locations can Inhibit vaporization at higher temperature locations. 

This distributes the loss of material over a wider area. The lifetime of a 
system would thus be longer than that predicted from simple Isothermal kinetic 
equations, perhaps by a factor of 1000. Some materials might thus be consid- 
ered for use In a system such as the reslstojet that might otherwise be ex- 
cluded from consideration. 


CONCLUDING REMARKS 

Throughout this report we have calculated rates In terms of g/cm^ sec 
and cm/sec. We have also expressed rates In terms of lifetime. Another use- 
ful way was In terms of Z, a relative rate with respect to the normal rate of 
vaporization. Actual vaporization rates are highly dependent on the proper- 
ties of the material. In contrast, Z values are not grossly dependent on the 
properties of the vaporizing material. Note that In general the Z values 
for Re throughout this report do not differ significantly from the respective 
values of Z for TaC. See particularly figure 19. Thus Z values calcu- 
lated for one material can serve as approximations for many materials. Once 
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the normal vaporization of a perspective material (e.g., zirconium dioxide) 
for an application Is calculated, its approximate rate under other conditions 
may be obtained by multiplying by the Z (for Re or TaC) under those condi- 
tions. The laborious calculations of Z from the equations given In this 
report are not required; an approximate value Is generally sufficient. 

APPENOIX 

Glossary 

a 1.01 3x1 0 6 /sqr (2irRT) 

(av) average 

b 1/RT 

(bl) boundary layer 

(cs) chemically formed species 

(ef) effective 

(h) cross section of hole through tube 

(l) phase 

(j) component 

(k) kinetic step 

k.j,k 2 , kinetic constants 

k 3 

(m) matrix gas 

n running subscript In Infinite series 

n(j) mole fraction of component 

(nv) normal vaporization 

(p) pressure 

(pc) particular condition 

(s) Inside surface of tube 

(ss) slowest step 

(v) vapor 

(l) first kinetic step 

(2) second kinetic step 

2 

A area, cm 

A' term giving persistence of molecular direction after collision 

B(n) term required In Infinite series; defined In equation (65) 

C carbon monomer 

C' activity of vaporizing component 

C"(s) concentratlqn of species diffusing In solid, g/cm^ 
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C"(v) 


D 

D(g) 

D(s) 

E 

E' 

F 

F' 

F(r) 

G 

G' 

K 

L 

M 

H' 

Model 

N 

N' 

N" 

P 

P' 

P(v) 

P(Y) 

Q' 

Q" 

R 

R' 

S 

S' 

S" 

T 

T' 

T" 

V 


concentration of species diffusing In gas, g/cm 3 
carbon dimer 
carbon trlmer 

o 

density, g/cm 

gaseous diffusion coefficient, cm 2 /sec 
solid state diffusion coefficient, cm^/sec 
see equation (66) 

Initial weight fraction of vaporizing species 

p 

flux of gas, g/cm sec 

3 

gaseous flow rate, cm /sec 

fraction of vaporizing component remaining In solid 

vaporization rate, cm/sec 

temperature gradient, deg/cm 

equilibrium constant 

length of tube 

molecular weight of vapor species, g/mole 
molecular weight of species lost from solid, g/mole 
see end of glossary 
Avogadro number, molecules/mole 
vaporization rate, moles/cm 2 sec 

3 

molecules/cm 
gas pressure, atm 

vapor pressure of pure component, atm 
vapor pressure, atm 

pressure of gas as function of Y, atm 

fraction of molecules that adhere to a surface upon collision 
molecular diameter, cm 
gas constant 

Inside radius of tube, cm 
Initial thickness of tube, cm 
boundary layer thickness, cm 
coating thickness, cm 
absolute temperature, K 
time at temperature, sec 
half-life, sec 
velocity of gas molecules 
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w 

X 

Y,(Y) 

Z 

Z' 

Model 


0 

1 

2 

XI 

X'l 

X"1 

X2 

Cl 

C2 

El 

E'l 

E2 

E ' 2 
FI 

F2 


2 

vaporization rate, g/cm sec 

distance out from the surface, cm 

distance from tube entrance, cm 

relative vaporization rate, W(pc)/W(nv) 

relative vaporization rate based on lifetime, T(nv)/T 

all numbered models are concerned with Incongruent vaporization; the 
rate of solid state diffusion (SSD) Is always a rate-controlling 
step. The letter and number following the word "Model" Indicates 
the other rate-controlling steps Involved; their meanings are 
listed below. 

only SSD 

SSD plus normal vaporization proportional to cube of surface 
concentration 

SSD plus normal vaporization proportional to first power of surface 
concentration 

SSD plus chemical vaporization proportional to square of surface 
concentration 

Model XI plus boundary layer diffusion 
Model X'l plus temperature gradient 

SSD plus chemical vaporization proportional to first power of 
surface concentration 

SSD plus surface coating effect; dependence on cubes of surface 
concentration 

SSD plus surface coating effect; dependence on first power of 
surface concentration 

SSD plus evaluated enclosure effect; dependence on cube of surface 
concentration 

Model El Including dependence on Y 

SSD plus evaluated enclosure effect; dependence on first power of 
surface concentration 
Model E2 Including dependence on Y 

SSD plus flowing gas effect; dependence on cube of surface concern- 
tratlon 

SSD plus flowing gas effect; dependence on first power of surface 
concentration 
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i REM TABLE 1(a) MODEL 0 

- 01-13. 9 8 REM dens i t y 

4 SI— . 9 3. fi R E M t h i c k n e s s o f rn a 
terial 

5 W 1 =. 000 "43 a R2=.8279 a W=W 1 
^(!tR2) a REM eive rate of v 
apo r i za t i on i n a/cfii2-sec 

6 M2=19£ 36 a M3 =12. 01 8 REM m 
o I e c u 1 a r weights 

7 02= . 6000001 1 3 a REM solid st 
ate diffusion coefficient. 

8 F 0 = 2 * 0 2 $ D 1 *: M 1 •••' S 1 / M 2 a REM F-r 
oduct of parameters needed f 
or infinite s e r i es 

10 71 = 100000000 1 $ REM initial t 
ime used for calculation 

11 F‘1=P1*PI 

1 99 S9=- ■: . 2 5 f. D 2 / < S 1 *S 1 > ) 

11.0 S 8=0 a s 2 = 0 a S3=e a N 0 =-i a 
REM start of calculation to 
r t i m e T 1 
Hi F2 = S9'.f:j 1 

113 N 8 = 1 a T9- 8 a REM N8 is the 
number of terms. added each t 
ime; T9 is the value of the 
last ter m 

1 1 y N 0 = N y + N y a R E M count ins t e r m 
s in t h e full s e r i e s 
1 — 1 REM line 122 bx passes soriie s 
teps when more t h a n one t e f rn 
is added t o the series at o 
n e t i r,i e 

122 IF N3<>1 THEN 187 

123 N2= I NT < . 5*N0> 

1 24 N4=<2*N2+1 >*<2*N2+1 * p i 
131 T 2 = F 0 : E X P (. F 2 % N 4 ') 8 REM term 
ip series for Z with halt nu 

rn b 6’ r o 1 t e r rn s 

18 7 M 3 = < 2 :* N 8 + 1 ) t ( 2 * h 8 + 1 > % p 1 
188 T 3 = 8 s H 3 * E X P < F 2 % N 3 > 8 REM ter 
m in series t o r Fir) w i t h t o 
tai terms 

198 T0=F0*EXP<F2*N3> 8 REM term 
i n full series to r 2. 

191 IF T0<1 . E-499 THEN GOTO 3l5 

192 S0=S0 + N8*T0 8 S2 = S2+ f=3 N R * T 2 

8 S3=S3+T3*N3 8 REM add i ns*" 
the ter m s to the series 

193 IF A3S<T9-T0> < . 01.*T0 THEN GO 
SUB 60 1 

194 REM line 193 tells pro-gram t 

0 a v e r a •? © if adoacient t © r m s 
•are close.; line 195 assures 
10 terms 

195 T9 = T 0 8 IF N 8 < 1 0 THEN 119 

i 9 6 D 4 = R B S v < S 8 - S 2 > / < $ © + 1 . E - 4 R 9 > t. 
1.97 IF 04 <.80G1 THEN 315 
t 9 REM lines 1 9 b — 1 9 7 tell f- r o s r 
a m t o s t o F' a d d l n -g ter m s if n 
ait series almost equals fui 

1 series 


199 

300 

315 

319 

320 

321 
526 
558 


559 

568 

600 

681 

604 


605 

700 


GOTO i 1 9 

R EM lines 315-526 p r i n t c a 1 c 
uated '/a lues 
PRINT i, T'= ,, ;T! 

P R I N T " SDR T 1 = " .; S 0 R ( T 1 '> 

PRINT "Z > REL Rfl= u . S0/W 
PRINT "F< r > .. BULK C0N = ";S3 
PRINT 

T 1 = T 1 x S Q R ( 3 0 } 8 REM detgrwi n 
es next time value for which 
to calculate the infinite s 
e r i es 

REM line 568 tells program w 
hen to end 

IF Tl<. 008001 THEN 780 
GOTO 318 

IF N8 < 1 08 THEN RETURN 
N0=N9-4.5*N8 8 N8=i0*H8 8 RE 
M when difference in terms i 
s small; 18 terms are aver as 
ed 

RETURN 

END 
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1 REM TRBLE 1(b) MODEL l 

2 ‘r’ 1 = 1 2 REM value of Z.; 028 

for Model Bl.; 0025 for El, 

. 00267 for Pi.: .00000143 for 

Cl 

3 REM change the value of Vi i 
n line 2 for the various mod 
els 

4 HI = 10 i? REM specifies 10 sli 
ces 

5 DIM C0< 1 1 > , ftl < 12> 

6 REM CO refers to content of 
slice.; ftl refers to amount t 
ransfered 

? Rl=1.8 2 REM ratio ot size o 
f slice HI to Nl + 1 

3 M2=192 .96 2 Ml=12.0i 2 REM m 
oiecular wts. 

9 01=13.9 2 REM density 

10 Hi =.009748 2 REM rate of vap 
orizaton under normal condit 
ions in g / c m 2 - s e c 

11 02 = . 0000001 1 l 2 REM solid st 
ate di f f usi on coe f f i c i ent cm 
2 s 0 c 

1 £ S0= .01 G REM t h i c kness of ma 
t e r i a 1 

13 T5=0 G REM setting ot time t 

0 zero 

14 ft 1 < 1 ) =0 G fll < 2>=0 a fll(Hl+ 2 ) 
=0 a REM setting original tr 
an$f er v slues 

15 T6= . 0@ 8 0 00601 & REM initial 
setting ot the time data sho 
u 1 d be p r i n ted 

lb C 0 <. 2 > ~ 1 a REM set tins a cone 
en t ra t i on 

17 ft 2 = 8 a REM sets the ini t ia 1 
amount lost at zero 

12 REM line 28 sets t h e vapor iz 
at to r, rates o t 4 v apor s p e c i 
e s 

20 G 1 =W 1 *M2 -'D 1 y M 1 2 G2 = G i * . 0279 

2 G3=G 1*574 2 G4=2.8*Gl 

21 G 3 = 0 2 G4= 0 2 REM delete 1 i n 
e tor Models XI, X‘l, and X" 

1 

26 REM lines 28 to 31 set the i 
n x 1 1 a 1 c one of all slices a 
nd transfers; 

27 REM and deternines the thin ten 
ess ot slice S 1 > at surface 

28 81 =000001 2 S 2 =S 1 2 S3=S2 

29 FOR N2=3 TO Nl+1 

30 C 0 ( N 2 > = 1 2 fll<N2> = 0 2 82=82* 
RJ 2 83=83+82 

31 NEXT M2 2 81=81*80/' S3 

34 T3= . 5*Sl*Si*<. 1+R1 >*RlxD2/2 2 
REM set a maximum time incr 
emen t 


85 


36 


37 

38 

39 
42 


43 

44 


45 

46 

47 
52 


53 

55 

56 


D f 
60 

6 1 
62 


63 

64 
166 

167 

168 

169 
1 70 

174 

175 
180 

184 

185 


188 


189 

190 


200 


T0= . 5*T6 2 REM lines 35 to 4 
9 select the time increment 

< T 0 > 

T4= . 5*8 1 y < C0 < 2 > *C0 <. 2 > *G 1 + C0 C 
2 > *G3+G2+G4 > 

IF T 4 > T 3 THEN T4=T3 
IF T4 <T 0 THEN T0=T4 
IF T5+T0M6 THEN T0=T6-T5 
REM line. 44 sets the amount 
of material lost at surface 
during (70) time increment 
ft 1 = < C0 < 2 > *C0 < 2 > *G 1 +C0 < 2 > *G3+ 
G2+G4>*T0*C:0':'2>*Y1 
S2=S1 2 REM lines 44 to 4?’ d 
et ermine the amount of trans 
fer between slices in time < 
T W > 

FOR H2=2 TO N1 

fil <. H2 > =2*D2*T0* < C0 (. N2+ 1 >-C0< 
N2 > ) y S2y C 1 +R 1 > 

S2=S2*R 1 2 NEXT N2 
REM lines 53 to 5? determine 
the new concentrations in e 
ach slice after time increme 
nt < T 0 > 

82=81 2 03=0 1 *T0 
FOR N2=2 TO Nl +1 
C0c:N2>=C0<N2> + <ftl <N2,t-ftl <;N 2 - 
1 > > y S 2 

S2=S2*R 1 2 NEXT N2 
REM line 61 corrects slice 1 
for amount lost via vapor iz 
at ion in time increment < T 0 > 
C9 ». 2 > =C0 < 2 > -ft 1 y'S 1 
2= V 1 *C 0 < 2 > * < C0 < 2 > *C0 2 > *G 1 +G 
2+G4+G3*C©<2> >y (Gl+G2> 2 REM 
calculates Z for each time 
CT5) 

ft2=ft2+fll 2 REM calculates to 
tal amount lost in time <T5> 
T5=T5+T0 2 REM calculates to 
tal time 

REM line 167 selects print in 

g or not printing 

IF T5 <T 6 THEN 35 

REM lines 169 to 180 print o 

u t calculat ed values 

PRINT " T ’ , TIME=" T5 

PRINT “ SQR T'=" ;SQR<T5> 

PRINT "Z,REL Rft=“.:Z 
PRINT "C 1 , SUR CON=*'.iC0<2> 

PR I NT “ F v R > , BU CON* " 1 -fiS/SO 
PRINT 

T 6 =T 6 *SQR<SQRt. 1 ©> > 2 REM set 
s new time for printing 
REM line 189 decides whether 
to end program 
IF T 6 > 100000000 THEN 200 
GOTO 35 2 REN add another ti 
me increment <. T 0 > 

END 
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1 REM TABLE Itc> MODEL 2 170 

2 VI si »? REM value of Z; .023 17 1 

for Model 62: 025 for E2; . 

00267 for F2 .00000143 for 
C2 

3 REM change the value of Y1 i 172 

n line 2 for the various mod 173 

els 174 

4 P2=. 0279+2 .3+574 G REM this 175 

term add effect of chemical 
contribution to normal vapor 176 

iration 177 

5 R2= . 0279 6 REM this term onl 
y adds effect of monomer vaF 

or iza t ion; delete this line 185 

for V.2 186 

8 M2= 1 82 . 36 G M 1 = 1 2 . 0 1 i» REM m 
circular wts . 

9 D 1 = 1 3 8 it REM density 138 

I 0 W 1 = . 000 74 8 G REM rate of V3F 

orizit ion un de r nor m a 1 c ond i 190 

t i ons : •? cm2 -sec 

II D2= . 0009001 1 1 G REM solid st 

ate d i f ♦ us i on coe f f i c i en t 191 

12 3 1 = . 0 1 G REM thickness of rna 192 

t e r i a 1 

13 W=Wi*<l+R2> G Z 1 = v 1 +P2 > / 1 . 02 193 

79 i? REM converts rate into 

terms ot total vaporization 
r a t e 

14 E = V 1 k 3 1 * W * M 2 • D 1 / 0 2 s M 1 I* E1=E 194 

4E it REM product ot paramete 195 

rs needed for infinite serie 

s calculation 

15 B 1 = 0 »2 63=0 it H3 = 0 G §9=02' < 196 

S 1 tS 1 > i © REM ass i ?n initial 197 

. values to variables 198 

16 T0= 1000O0000 G REM initial t 
i me 

17 3=0 it $2=0 it S3=0 S N0=-1 G 

N8 = 1 i? T5 = U 8 T=T0*S9 G REM 19? 

initial parameters for each 

time parameter 200 

18 REM NS is the number of term 300 

s added to series at one t i ra 

e 310 

120 N0=N0+NS it REM count of term 315 

s 316 

122-82 = 1 it REM first evaluation 319 

for B 1 328 

130 REM lines 14Q to 151 evaluat 321 

e B 1 Fa raw? ter for infinite 526 

series 550 

140 81= 9*82+ 1*81 

150 B2=RrHCE/Bl >+N0:*:PI 

151 IF RBSCB1-B2.>> . 0000001*82 TH 

EN 140 559 

169 REM line 170 bypasses some s 

t eFs when the average of sev 560 

oral terms are added to inf. 599 

series 


600 

605 


649 

700 


IF N8 0 1 THEN 185 
REM lines 172 to 176 calcula 
te 84 for infinite series wi 
th only half the number of t 
e rms 

N2=INT< . S*H0> G 84=1 

83= . 9*84+ . 1*83 

B4=fl TN*: E'*B3 >+N2 fcPI 

IF flBS'%B3-B4>> .0800001*84 TH 

EN 173 

86=84*84 

T2=N8*EXPt-' B6*T>)/<B6+E1+E> 
•2 REM term for half series 
for Z 
85=62482 

T4=N8 *EXP < - < 65*T ) ;»/< B5+Et+E> 
it F:EN term for full series 
for Z 

T3=T4 G5 i? REM term tor full 
series tor* F < r > 

S = S + T 4 G S2=S2+ 5*T2 8 33*93 
+T3 8 REM addins terms to fo 
rrn series 

IF T4 •• 1 . E-4 0H THEN 31E 
IF flBS<T4-T5^ < . 01*T4 THEN GO 
SUB 660 

T 5 = T 4 it REM lines 192 & 193 
dec x de whe t he r t e r n» i s sma 1 1 
enough to take the average 
of terms 

IF N 0 < 1 0 THEN 120 
REM line 194 makes sure that 
the fist 10 terms are not a 
ve rased 

03=ABS< <8-92,’/<S+l . E-50) > 

IF D3 < . 000 1 THEN 310 
REM lines 196-197 tells pros 
ram to stop adding terms if 
half series almost e^ual ful 
} series 

REM line 200 says to go back 
a nd calculate nex t t e rms 
GOTO 120 

REM lines 310-526 say to pri 
nt the calculated values 
C 1 =2*S*E G F9=2*S3*E1 
PRINT M T' = " .5 T0 
PRINT "SQR T 1 = ,# :• SQR<T0> 

PRINT "Z.REL RATE=" s C1*Z1*Y1 
PRINT *‘C *. .• SUP C0N= H ; Cl 
PRINT "F< r > • BULK C0N= M ;F9 
PRINT 

T0=T0-SQF:< 10> ? rem de term in 
es next time value for which 
to calculate the infinite s 
er ies 

REM line 560 tell program wh 
en to end 

IF TOC. 00O001 THEN 700 
GOTO 17 © REM go back and st 
art calculation on next time 
parameter 

IF N0\ 1O0 THEN RETURN 
N0=N0-4.5*N8 G N8=M8*18 G RE 
M when difference in terms i 
s small. 19 terms are averag 
ed 

RETURN 

END 
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Temperature, 



Vaporization rate of rhenium, g/cm 2 sec 





10 ‘ ; 


g 10' 4 



10' 6 

II 


Coating thickness, S", cm 
Figure 5. - Effect of coating on vaporization. 



Figure 6. - Effect of solid state diffusion on the composition of TaC during vaporization [thickness 
0.01 cm; symbols on curves refer to models]. 




Relative vaporization rate, 


Re in flowing gas- 


TaC in flowing gas 
(normal vap) 





Average bulk concentration of carbon, F(r) 




l(f 3 10* 2 10° 10 2 10 4 10 6 10 8 


Time, sec 

Figure 13. - Combined effects of enclosure and diffusion due to in- 
congruent vaporization on carbon concentration in TaC[ numbers 
on curves areY values (cm) for Model E'2; thickness - a 01 cm; 
A(h)/A(s) - 0.00251 


Relative vaporization rate, 


m 


Time, sec 

Figure 14. - Dependence of vaporization of Ta upon combined effects 
of enclosure and diffusion assxiated with incongruent vaporiza- 
tion [numbers on curves are Y values (cm) for Model E’2; thickness 
-0.01 cm ; A(h)/A(s) - 0.00251 


Distance from entrance, cm 

Fl [M^ 1 | 5 E2i C ° mbined effectsofenc,osure and dif fusion due to incongruent vaporization on half life of TaC 



Figure 16. - Interaction of flowing gas in tube with bound- 
ary layer and coating effects. 



Value of parameter indicated on curve 



Figure 17. - Interaction of stationary gas in tube with 
boundary layer and coating effects. 



Figure 18. - Calculations for the vaporization ofTaC involving chemical reaction and boundary layer 
diffusion. 





Distance from entrance, cm 


Figure 19. - Vaporization rate profile for a gradient of 10° 
per centimeter. 



PM 


8 - 


Distance from entrance, cm 

Figure 21. - Effect of various parameters on the veporiza- 
tion profile of TaC with 10° per centimeter gradient 



Time, sec 


Figure 22. - Calculations for the vaporization of TaC involving temperatui 
gradient, chemical reaction, and boundary layer diffusion [G' - 10°/cm* 
S' - 0.005cm]. 
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